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Heat Storage in Alloy Transformations 


Introduction 

This project has explored the eutectic transformation 
of alloys containing plentiful metals as a means of storing 
heat at moderately high temperatures. In the course of this 
work, the theory of eutectic transformation was examined 
to find guidelines to the best material combinations to 
examine. The heats of transformation were measured 
calorimetrically , and the volume changes of expanding 
solid mixtures and homogeneous liquid solutions, and 
especially during the transformation between the two 
states at fixed temperature, were measured by changes in 
X-ray absorption. Heat flow models appropriate to storage 
in phase change materials have been developed along with 
efficient calculating procedures so that the relative 
importance of the problems associated with energy storage 
density, heat conduction and similar properties can be 
assessed in light of what is possible now. In the plan 
to use alloy eutectics as heat storage media the one step 
that is not completely resolved and demonstrated is 
containment in a manner that will retain the advantages 
of excellent thermal conductivity. It has been demonstrated 
that silicon carbide resists attack by the molten storage 
alloys. However, it has not been demonstrated that carbide 
films can be grown on oxidation resistant alloys and 
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retained through many melting and freezing cycles of 
the storage system. It has been proposed to extend the 
research work in this direction so that working storage 
systems can be designed and tested once containment is 
assured. 

This report is constructed with the cited 
conditions in mind. The first section is a summary 
statement intended to show what has been accomplished in 
preparing heat storage systems utilizing alloy eutectic 
transformations as the working core. The discussion is 
in fairly general terms, because the detailed support 
for the statements is given in a series of extended 
appendices. The appendices are reprints of pub 1 .shed 
papers or papers being prepared for publication. In 
several cases more extended discussion of the topics is 
available in Master's or Ph.D. theses accepted by the 
University of Delaware or which soon will be submitted 
to the faculty of the University of Delaware. 

Heat Storage in Eutectic Alloys 

Appendix A, "Heat Storage in Eutectic Alloys," by 


C. Ernest Birchenall and Alan F. Riechman, Met. Trans. A, 
llA , 1415-1420 (1980) , contains the thermodynamic equations 
governing eutectic transformations in binary systems. The 



relative costs of elements and their availability indicates 
a current preference for S, Al, Si, Zn, P, Na, Cu, Mg, 

Ca. Fe, although plentiful and cheap, is most likely to 
be needed in the solid containment systems. 

Metal or alloy heat storage by phase change 
appears to be feasible from 1219K to moderately low 
temperatures. Inorganic salts with comparable storage 
capabilities per unit mass or volume also are available, 
as shown in Table II, Appendix A, but the thermal 
conductivities are one to two orders of magnitude higher 
for the metallic storage materials than they are for the 
salts. Simpler and smaller heat exchangers should be 
possible with the metals than with the salts. 

The theory of heats of eutectic transformation shows 
that the important contributions come from heats of fusion 
of the participating elements, possibly modified by their 
combinations in intermetallic compounds, and from the 
differences in the entropies of mixing of the relatively 
well-ordered solid phases and the eutectic liquid. 

Alloys were chosen for study, or sought where they were 
not previously known, in order to optimize the heat storage 
density. For a current tabulation see Appendix B. Although 
it is not possible to state that the best alloys have been 
found, reasonable assurance can be given that the reported 


alloys are representative of good thermal storage behavior. 
They give the alloys a fair test in comparison with other 
storage materials. 

It is important to consider one other cost factor 
which appears to favor the alloys. Even after 
prolonged use, they remain refined metallic materials 
that have other valuable uses. They retain much of their 
initial value. As a specific example, the Al-Si eutectic 
alloy has been widely used in cast automotive parts, 
including engine blocks. Although many other alloys are 
not in precisely the compositions that would be used in 
other applications, they could be adjusted by simple 
additions. 

Volume Change During Transformation 

To maintain efficient heat transfer from the PCM 
containment vessel to the heat transfer fluid or to a 
heated surface such as the receiver cavity of a solar 
receiver the solid or liquid PCM must always be in direct 
contact. Large volume change during transformation will 
require design complexity to insure that direct conducting 
paths are maintained. The large volume changes reported 
for halide eutectics, on the order of 20 to 25 percent of 
the total volume, were sufficient incentive to measure the 


changes that occur for the metal eutectics. 


An X-ray absorption method that is readily used 
with ductile samples is described in Appendix C. Its 
application is illustrated by measurements on pure aluminum 
and the Al-Cu and Al-Si eutectics that agree well with 
previously published values. Thermal expansion 
coefficients for the liquid and solid phases or mixtures 
of phases also are in reasonable agreement where published 
values are available for comparison. 

Some refinements have been added to the method in 
Appendix D to simplify the selection of a useful X-ray 
wavelength and improve control of furnace temperature. 

The measurement of fractional volume change during 
eutectic transformation has been extended to two ternary 
systems, with work still progressing. Expansion coefficients 
have been measured again. 

The alloy eutectics show volume changes of about 
5 to 7 pet. for the aluminum alloys and near 3 pet. for 
a Cu-P-Zn eutectic. There is no reason to believe that 
these modest changes are atypical, so the operation of a 
PCM storage system should be simpler with these materials 
than with those that have much larger volume changes 
during transformation. 


Modeling Heat Flow in Cylindrical Heat Exchangers 
Containing Phase-Change Materials 

Appendix E describes the development of a 
numerical routine for calculating heat flow efficiently 
in a vessel of cylindrical symmetry containing a material 
capable of storing latent heat during a solid to liquid 
transformation, including the special case in which the 
solid is a mixture of phases. This paper was presented 
at the 3rd Miami International Conference on Alternative 
Energy Sources, December 15 to 18, 1980, and will appear 
in the proceedings of that conference to be issued by the 
Hemisphere Publishing Company. The utility of the procedure 
for dealing with metal and salt eutectic storage materials 
and the factors affecting their performances in simple 
heat exchangers is demonstrated in this work. 

An improved procedure is outlined in Appendix F 
that has all of the capabilities of the earlier procedure 
but is capable of more precise calculation. Results have 
been collected with this procedure right up to the 
preparation of the report, with additional runs to be 
made. The paper is being prepared for the AIAA Thermophysics 
Conference, Palo Alto, Cal., to be held in June, 1981. 

The calculations show the advantages for the highly 
conducting alloys in giving faster storage and retrieval 
of heat for fixed geometry as compared with those materials 
that have lower thermal conductivities. 


These limited applications have given assurance 
that the calculation techniques can be used in comparing 
alternative designs as well as competitive PCMs in 
establishing optimum system designs. Following demon- 
stration of a proven containment procedure, heat exchanger 
designs consistent with the containment method and storage, 
materials will be evaluated, leading to the choice of the 
optimum configurations for the various applications that 
seem reasonable for these combinations. 

Containment 

Recent studies that show the short term stability 
of silicon carbide in contact with molten storage alloys 
are described in Appendix G. Alloys rich in aluminum 
or copper, containing significant concentrations of 
silicon, magnesium and zinc have shown little tendency 
to attack silicon carbide, but an appreciable rate of 
attack on Type 304 stainless steel (19Cr, 9Ni) . Although 
the exact properties of the containment materials that are 
desirable will be determined by the geometrical con- 
figuration selected, it seems likely that the structure 
will have to be fabricated from an alloy that combines 
good high temperature strength, good thermal conductivity, 
and reasonably good low temperature ductility. In addition, 
resistance to chemical attack by the heat storage materials 
must be provided without impairing these other properties. 


One reasonable procedure would be to construct the 
system from a stainless steel or other oxidation resistant 
alloy, then coat surfaces that will be exposed to the heat 
storage alloy with silicon carbide. The coating would 
cause little loss in heat transfer and would provide a 
small increase in high temperature strength. If the 
surfaces exposed to air were also coated, somewhat lower 
grade alloys might be used, because silicon carbide should 
provide good oxidation resistance. 

It is important to emphasize that this item is the 
last bit of undemonstrated technology needed to embark on 
a serious design. In particular, the stability of a 
carbide coated steel to many thermal storage and release 
cycles should be verified on a laboratory scale before the 
construction of a working system is undertaken. Several 
substrate alloys, differing in their thermal expansion 
properties, should be coated and tested before a final 
choice of containment material is made. 

A more complete summary of Appendices E and F 
prepared by Michael Labdon follows. 


INTRODUCTION 


The thermal storage units depicted in Figures 
1 , 2, and 3 represent conceptual designs for use in the 
collection and transmission of energy concentrated from 
parabolic-dish solar receivers. Due to the high energy 
flux imposed on the inner cylinder, the ongoing 
investigation dealt with the contrasting of candidate 
phase-change materials (PCM's) for conducting the heat 
energy, as well as finding an appropriate material for 
PCM containment. 

Due to the array of potential PCM candidates and 
the complicated geometries associated with the storage 
structure, an efficient and versatile numerical routine 
was required for analyzing their thermal behavior. 


PROCEDURE 

Method 

An extensive literature search was initiated to 
determine if a solution method already existed which 
could be modified, or if an original routine needed to 
be developed. The result of the search was to utilize 
a complex method originally posited by Forster (1) , 
Murray (2) and Lazaridis (3), which was subsequently 
demonstrated as workable by Farouk (4) , and finally 
clarified by Harding (5) and Labdon (6) as to its 
inherent simplicity and efficiency. 
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The versatility of the routine is reflected in the 
fact that time and spatial dependent boundary conditions, 
temperature-dependent properties and initial temperature 
distributions may be included, which represents the more 
realistic condition encountered with solar energy and 
storage systems. Outputs include prediction of melt-line 
locations, temperature histories and charging rates. 

The method models heat conduction within an annulus of 
phase-change material subjected to various boundary 
conditions, by overcoming the inherent difficulty of 
determining the temperature distributions around the 
melt-line (hence, its velocity) which is non-linear in 
nature. Essentially, finite-difference equations monitor 
the heat transfer in the regions away from the interface, 
while Taylor-series expansions of temperature are made 
near the fusion front so as to describe the temperature 
fields and thus couple the thermal behavior between the 
solid and liquid regions. 

The mathematical derivations, assumptions and 
solution process of the numerical routine appear in two 
publications, Appendix E and Appendix F, where the latter 
models both the radial and axial heat conduction within 
the material. The various thermal properties of the 
candidate PCM's, and the assumptions used in facilitating 
the generation of results are listed on page 15, Appendix E. 


Analysis of Storage Units 

The solar receiver/storage units dimensioned 
in Figure 1-3 do not represent an attempt in design 
optimization , but merely provide a parity for the contrasting 
of any desirable phase-change materials, which in this 
case are metal alloys and salts. 

Figure 1 was subjected to the simultaneous 
application of heat flux and convection on the inner and 
outer surfaces respectively, with the entire medium 
containing a 10°C initial superheat. In this manner, an 
evaluation of the PCM's can be made while the storage unit 
is in the actual working mode. 

Figure 2 differs from that of Figure 1 by the 
cylindrical array of internal tubes which contain the 
heat-convecting medium. Only the discharging case will 
be modeled by assuming the entire medium is initially in 
a constant molten state. Six equally-spaced, 4 cm. diameter 
tubes are situated halfway between the inner and outer 
surfaces. Termination of the discharging case occurs 
when the radially-expanding melt-lines intersect each 
other and the unit's surfaces. The small amount of material 
left over is taken into account by forming an extra concentric 

ring around each tube. 


The third storage unit, depicted in Figure 3, 
represents a more complicated configuration presently 
outside the scope of this work, but it is reasonable to 
expect that the contrasting performance between the salt 
and metal alloys would be similarly evident in this 
geometry as well. 


RESULTS 

Figures 4 through 9 represent output from the 
routine for the first solar receiver/storage unit 
(Figure 1). Figures 4, 5, 6 represent melt-line locations 
within the receiver for the salt, Al-Si, and Mg-Si data 
respectively. The melt-line contours actually represent 
a three-dimensional surface similar to a cone or trumpet- 
end. The more curved lines of the alloys represent the 
increased heat transfer within the upstream nodes due to 
the higher thermal conductivity. The smaller times 
required for melt-line penetration are also indicative of 
the alloys. Interface locations are shown only for the 
outer half of the unit due to the large computer process 
times required for output generation, and due to the fact 
that a steady-state operation is approached where the heat 
flux imposed on the inner cylinder equals the heat 
convected out, thus yielding a melt-line velocity nearing 


zero. 


Figures 7 and 8 describe the heat flux boundary surface 
temperatures for equal time increments. Both show the 
steady-state temperature distribution being approached, 
as well as the much higher surface temperatures 
associated with the use of the salt PCM; attributable 
again to the low thermal conductivity. Figure 8 indicates 
the effect of increased melt-line penetration on the 
"cooling" or surface temperature reduction of the unit, 
typical of metal alloys. 

A measure of the effectiveness of this type of heat 
exchanger is gleaned from the charging rates shown in 
Figure 9. The area under the curves represents the 
amount of energy transferred within a certain amount of 
time, and as can be seen, the metal alloys are more 
efficient in storing and transferring the collected energy. 

As a further display of the versatility of this 
routine, the discharging rate of Figure 2 was determined 
by the amount of energy transferred from the six tubes, and 
contrasted to that obtained from the discharging case of 
Figure 1 (when the heat flux was terminated) , the results 
of which appear in Figure 10. This indicates the role 
the method plays in design optimization of the receiver/ 
storage units, by showing the increased discharging rates 
associated with installing the convection tubes within 
the PCM, as opposed to a sheet surrounding the entire 
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ABSTRACT 


Eutectic alloys and intermetallic compounds of the 
common elements, Al, Ca, Cu, Mg, P, Si, Zn, may be useful 
for high temperature heat storage. A number of binary and 
ternary eutectic alloys have already been studied and their 
heats of fusion measured. In this work heats of fusion, and 
in come cases heat capacities, of new multicomponent 
eutectics and intermetallic phases are reported. Some of 
the alloys investigated are competitive with molten salts 
for heat storage density at high temperatures. 

A method to determine unknown eutectic compositions 
was developed, based on combined results of differential 
thermal analysis with metallography and microprobe analysis. 
A new differential thermal analyzer, which gave satisfactory 
accuracy for qualitative measurements, was built to handle 
the relatively large samples that are required for this 
procedure. The method allows determination of eutectic 
composition in no more than three steps. 

The heats of fusion of the resulting alloys were 
measured using commercial calorimeters, a differential 
thermal analyzer and a differential scanning calorimeter. 

The possibility of measuring heat capacities by 
differential thermal analysis was also examined. 


INTRODUCTION 


There are various ways in which heat storage can 
reduce costs and make more efficient use of power 
generating equipment. A thermal storage unit would 
eliminate the difference between peak and off-peak hours, 
storing heat when demand is less than the generating 
capacity and returning it as needed. The generator could 
be designed for the average, instead of the maximum, power 
requirement. A central storage system could thus reduce 
the cost of building and operating a power generating plant 

An efficient way to store heat is to use phase 
transformations that take place at a definite temperature. 
The storage density depends upon the latent heat of 
transformation, which is the entropy change times the 
transformation temperature. This entropy change is small 
for solid-solid transformations but much larger in melting. 
For solid- gas transformat ions, advantages of the still 
larger entropy change are offset by the practical problems 
in handling the large volume of gases. 

The first materials considered for high temperature 
heat storage were inorganic salts. Fluoride eutectics 
were studied by Schroeder (1). Lithium compounds and 
eutectics were characterized by Tye et al. (2) . Further 
study in molten salts was carried out by Ferrara et al. (3) 


by Haru et al. (4), and by Clayton (5). Some of the 
problems found in molten salts were: 

i) They have low thermal conductivity, leading 
to large irreversible losses and requiring 
a very large surface to volume ratio. 

ii) They are highly corrosive and must be dried 
very carefully, even if expensive, corrosion- 
resistant containers are used. 

iii) They have large volume changes during melting, 
around 20%. 

iv) Lithium and fluoride compounds, those with the 
best heat storage densities, are relatively 
expensive materials. 

Metals diminish these problems, and therefore might 
be strongly competitive even if they have somewhat lower 
storage densities. Materials were selected in the range 
of 300 to 1000°C, on the basis of cost and availability, 
as well as storage density (6) . The heats of fusion of 
many binary and ternary eutectics were calculated 
approximately (7) , and the most promising of these 
were measured (8) . 

In choosing particular eutectics, calculations were 
done for the heats of fusion using two different ideal 
solution models as well as a regular solution model. 


These calculations are discussed below. The results 
are only approximate, but it was possible to identify 
qualitatively the contributions to the heat of fusion 
in order of importance. 

i) Heats of fusion of the constituent elements, 

ii) Heat of formation of intermetallic compounds 
present. 

iii) Difference in entropies of mixing between 
solid and liquid phases. 

Analysis of the results obtained by Riechman shows that 
alloys have the best observed heat storage densities in 
the range from 700°K to 850°K. Better alloys are sought 
in the region 900°K to 1000°K and below 700°K by further 
all-tying. Therefore, it is necessary to determine the 
conyruently melting compositions in multicomponent systems 
and evaluate the heats of fusion. 

Ternary phase diagrams are not generally known, 
and few thermal data are available. In some cases 
the formation of ternary intermetallic compounds is reported. 
With this lack of information any theoretical approach to 
the problem is a crude approximation. Nevertheless, such 
an approach can be the useful starting point of a trial- 
and-error method for experimental determination. 


Although it is clear that alloys will be used primarily 
as latent heat storage systems and only incidentally for 
sensible storage, it is important to know their heat 
capacities as well. Complete information is needed to 
analyze the possible heat storage systems and develop 
practical design. 


CALCULATION OF EUTECTIC COMPOSITIONS AND TEMPERATURES 

For the ideal solution model all calculations are 
simple, and no thermal data are needed except heats of 
fusion of the components. For the regular solution case, 
the binary equations have been extended to ternary systems, 
an approximation that allows calculation of ternary free 
energies from binary data. This was first done by Darken 
(9) in 1950 and developed further by Toop (10) . The 
derivation is very well presented by Richardson (11) . 

An alternative but equivalent treatment is given by 
Olson and Toop (12). This approximation has been used 
successfully in the calculation of phase diagrams with no 
ternary interraetallic compounds, for example, the Ni-Cu-Cr 
system by Meijering (13) . Tne limitations of the method 
were analyzed by Olson and Toop (14) . In cases where 
compound formation is suspected, they concluded that the 
regular solution model does not yield useful results. 

In addition to these two equations, a polynomial 
approach has been used. This Guggenheim-type equation has 
been generalized by Hillert (15) . 

The three approaches were compared by Ansara et al. 

(16) for the determination of phase boundaries. In the 
systems Fe-Cr-V and Cd-Bi-Sb it was found that the results 
derived from these three equations did not differ much. In 


cases where compounds are formed or where there is a strong 
ordering in the liquid these approaches are not expected 
to be useful, especially the first two models that assume 
regular solution throughout the system. The advantage of 
these models is that they are based on binary data and no 
ternary measurements are needed. 

In this work both ideal and regular solution models 
were applied and the results compared. Clausius-Clapeyron 
equations, one for each component phase, and the condition 
that all the fractional compositions sum to unity yields 
n + 1 equations in n unknown compositions and the 
temperature. This set of non-linear equations was solved 
by computer fitting. The predicted eutectic temperatures 
for some alloy systems are shown in Table 1. 

A regular solution model applied on the basis of 
pure components will not give adequate results if 
intermetallic compounds are present in the system. For 
this reason the same approximation was used taking the 
constituent solid phases as the components of the system. 
Thermal data for the calculations were taken from 
Hultgren et al. (17), except in the Mg-Si case (18), 
where activities are always given with the pure elements 
as standard states. These reference states were changed 
to be the pure intermediate phases. The eutectic points 
predicted by this model are also shown in Table 3 . 




CALCULATION OF THE HEATS OF FUSION 
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For eutectic melting the ideal solution model can 
be applied to the prediction of heats of fusion using either 
the participating elements or the constituent phases as 
components. The two approaches yield different results. 

Both assume complete ordering in the solid with different 
ordering in the liquid. Random mixing of elements is 
likely to be better for disordered liquids and random 
mixing of phases for partially ordered liquids. 

Regular solution equations for eutectics (19) 
permit calculation of the latent heat either from the 
entropy change or from the enthalpy change , assuming 
random mixing in the liquid. If the solutions are truly 
regular and the data are correct, both equations should 
be equivalent. 

Table 2 is a comparison between the described 
theoretical methods and experimental results. It 
appears that the ideal mixing calculation from the parti- 
cipating phases often is the best choice for accuracy and 
simplicity. The calculations by all methods tend to 
overestimate the heat of fusion, and they should be taken 
as an upper bound to the true value. 
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EXPERIMENTAL DETERMINATION OF EUTECTIC COMPOSITIONS 

An experimental method for finding and characterizing 
eutectics was needed that would work with no information at 
all on the system but would be able to use known information 
to converge faster. The combined results of DTA, electron 
microprobe and metallography provide such a trial-and-error 
method. Any theoretical approximation or empirical data 
can be used to improve the initial guess of eutectic compo- 
sition. 

A sample of an arbitrary cr guessed starting composition 
is prepared, and the melting behavior is determined by DTA. 
Metallographic inspection of the alloy is done optically 
and in a scanning electron microscope (SEM) . The observed 
phases are correlated with the DTA results. The sample is 
searched for a small eutectic region, and the composition of 
any small zone suspected of being a eutectic mixture is 
analyzed quantitatively by SEM with an energy dispersive 
analyzer for the emitted x-rays (EDAX) . If some doubt 
remains about the number of identity of the phases present, 
x-ray diffraction measurements are made. 

If there is a eutectic region with a melting tempera- 
ture in a desirable range, a new sample is prepared with 
the composition obtained by EDAX on the previous sample. 

The procedure is repeated until an alloy that is completely 
eutectic is obtained. When an intermetallic compound is 


found in the process, it is also studied, especially in 
higher order systems where the probability of finding a 
desirable eutectic is lower. 

In order to be sure that the observed DTA peaks corres- 
pond to the observed phases and the compositions measured, 
it is necessary to work with the three techniques on the 
same sample. In this case inhomogeneities in the alloy are 
not a serious problem. 

k 

Commercial calorimeters v :k with very small samples in 
order to minimize the errors in quantitative measurements. 

A DTA apparatus that will give good qualitative measure- 
ments with a large sample is diagrammed in Fig. 1. Temper- 
ature was varied under control or a programmer. The 
differential temperature was amplified and plotted directly 
against sample temperature on an XY recorder. 



PREPARATION OF ALLOYS 


The first alloys were made in a resistance furnace 
following the procedures of Riechman (3,19). However, 
when the starting materials differed greatly in melting 
temperatures, vapor losses could not be avoided, even with 
a collar on the stirrer shaft closely fitted to the crucible. 

Owing to the repeated failure to prepare the 
alloys in which the starting materials have large 
differences in melting temperatures, an induction furnace 
was assembled. The relatively low frequency (50 kHz) of 
the mercury spark gap resonator provided vigorous internal 
stirring to the melt. 

The starting materials were introduced into alumina 
crucibles that were then sealed by an alumina disc of 
the same diameter cemented with Sauereisen cement No. 78. 

In the cases where most of the starting materials were semi- 
conductors, graphite crucibles were used, sealing them in 
the same way. The crucible was place;' in a mullite tube, 
with a thermocouple touching it from the bottom and argon 
flowing in from the top. This setup was entirely successful 
in preventing vapor losses and providing good stirring to 


the melt. 


MEASUREMENT OF HEATS OF FUSION AND 
HEAT CAPACITIES 

Two commercial instruments were used for the 
measurement of heats of fusion, the Perkin Elmer DSC-II 
and the Du Pont 990 DTA. The former gave better precision, 
about 3%, and was used for alloys with melting points 
be. w 700°C. The Du Pont 990 DTA with its high temperature 
cell gave measurements with 5% accuracy for temperatures 
above 700°C. 

For alloys with melting temperatures below 700°C 
the measurement of heat capacities can be carried out by 
DSC with a 10% accuracy. Manufacturers of commercial 
DTA instruments for measurements at higher temperatures 
do not claim sufficient baseline control for heat 
capacity determination. However, the Du Pont 990 DTA was 
tested for measuring heat capacities because of its 
improved electronic control of the baseline. Twelve 
baselines were analyzed, and it was found that they were 
straight within 10%. To check reproducibility several 
runs were made with the same sample pan that was opened 
and repositioned in the furnace between runs. 

After two successful calibrations with sapphire, 
measurements on Al and Ag gave results shown in Table 3. 
For Al the agreement with tabulated data is very good. 

For Ag agreement is not as good, because its heat capacity 


Table 3 


i: 

i. 

i: 

j Heat Capacity of Al and Ag 

Measured with DTA 1200°C Cell 


T(C°) 

Cp(meas) 

Cp (tabulated) 


(cal/g) 

(cal/g) 

Aluminum 

550 

• 

ro 

ON 

.28 

580 

.26 

.28 

600 

.26 

.28 

630 

.27 

.29 

650 

.29 

.29 

(liq) 710 

.36 



Silver 


800 

.12 

.07 

850 

.10 

.07 

900 

.15 

.07 

930 

.25 

.07 

(liq) 1000 

. 5 

- - 










taken at higher temperatures, is much lower. At lower 
temperatures the agreement improves. 

Measurements were also attempted for the Mg-Si, 
Cu-Si and Mg-Cu-Si eutectics. At high temperatures, the 
values obtained are too low. In some cases the solid and 
liquid baselines crossed, owing to the evaporation of 
the volatile elements. 

It appears that the specific heat values from DTA 
for non-volatile alloys are reliable within about 20%, 
but the values for volatile alloys cannot be trusted. 


RESULTS 


New Eutectics and Intermetal lie Compounds 


A search for new eutectics and intermetallic 
compounds was conducted in several ternary systems and 
two quaternaries. The search method converges in three 
steps, in two cases to negative results. The diagram 
in Table 4 summarizes the results step by step. 

X-ray data on the new compounds found are reported 
separately (20) , together with several interesting 
microstructures . 


Measurement of Heats of Transformation 


Heats of eutectic transformation were measured in 
the Perkin Elmer DCS II when the melting temperature was 
lower than 700 °C. For higher temperatures the Du Pont 
DTA model 990 was used with its high temperature cell. 


Several known eutectics were characterized. 

Two of the previously measured eutectics (19) were 
repeated, Mg-Zn because the previous result was believed 
to be too low, and Mg-Si for further check. For Mg-Zn a 
higher value was obtained, and for Mg-Si excellent 
agreement was found with the previous measurement. The 
results for the known alloys are listed in Table 5, 
together with those newly determined. Table 5 also 
includes the measurements of heat capacities for the 
alloys that are interesting as heat storage materials. 
The room temperature densities used in calculating 


able 4. 


Steps in the.- Search Cor New Eutectics and Intermetallic Compounds 




P-Cu-Si 



Alloy No. 1 j 

Microstructure Thermogram | 

M V i 7 Cj 16 + 

segregated 

Mg 2 Cu 

2 peaks 

ternary 
eutectic + 
segregated 
Mg 2 Cu 

1 peak 
453®, bump 
after t-.he 

MgSi cut. 4 
Zn 3% Mg 
eutectic 

1 broad 
peak 
000®C 

eutectic: 

8 brass, 
Cu.P, 

ternary comp. 

1 peak 
720®C 

1 hcmogeneou3 
phase 

1 peak 
84 0 ®C 

ternary 
compound + 
ternary 
eutectic 

2 peaks 

CaSi 2 4 Si 
4 CajMgSi^ 

2 peaks 

AljCu 4 Si 
4 Mg, Si 4 
quat. comp. 

4 peaks 

2 quat. comp. 
♦ quat. 
eutectic 

2 peaks 

Al-Mg-Zn eut. 
4 Kg 2 Si 4 

Si 

2 peaks 


Alloy No. 2 Alloy No. 3 

Mtorostructure Thermogram Microstructure Thermogram 


M ^6 S 17 Cu 16+ 2 P« ak8 

eutectic 


ternary 
eutectic: 
Mg 2 Ca, Mg. 

ternary comp. 


Mg-Si eut. 
+Zn 3% eut. 


ternary 

eutectic: 

Cu^Si, 

ternary comp, 
other phase? ■ 


CaSij + Si 

4- CajHgSij 


AljCu + Si 
+ M<J 2 Si 4 
quat comp. 


eutectic: 

2 quat .comp., 
Al, Other 
phase 


Al-Mg-Zn eut. 
4My 2 Si+Si 


ternary 

eutectic: 

Mg 6 Si 7 Cu 16' 

Mg 2 Si,Si 


1 peak 
445*C 


1 broad 
peak 
850*C 


Mg-Si eut. 
4 Zn 


2 peaks 


4 peaks 


1 peak 

52 0*<-' 


2 peaks 


Ilg 2 Si 4 si 

4 Ca 2 *MgSi 3 


1 peak 
770*C 



1 broad 
peak 
900°C 



1 peak 
865®C 











































Table 5. Summary of Heat of Fusion Measurements 
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DISCUSSION 


The search method exceeded expectations by 
converging in three steps in the most complicated 
case. SEM metallography was used in almost all 
cases. Even with electron microscopy it was not 
possible to distinguish easily all the phases present. 

In a number of systems the phases had to be recognized 
either by composition or X-ray diffraction measurements. 
Eutectic mixtures also were difficult to recognize; 
lamellar morphologies seldom occurred. In addition to 
metallography, comparison of the overall composition to 
that of the fine, clearly eutectic region was helpful. 
Size of the phases was not by itself an adequate 
criterion. Extreme cases were found. Mg-Cu-Ca had such 
a small eutectic spacing that X-ray diffraction studies 
were necessary to determine whether there was a eutectic 
or a single intermetallic compound. Several cases 
exhibited widely varied spacings in different regions 
of the sample, probably as a consequence of different 
freezing rates. 

Compositions by EDAX can be measured within 3% 
absolute error. The most serious limitation of the 
EDAX measurement is its spatial resolution. If a small 
particle of a phase to be analyzed is embedded in a 
different matrix, the penetration of the beam at the 





lowest possible voltage in some cases excites the matrix 
elements beneath the particle. Also it is impossible to 
get precise composition values for the different phases 
that constitute an eutectic if the spacing is small. 

The performance of the DTA for large samples 
was entirely adequate for qualitative measurements. 

A 30°C difference between the measured and actual 
transition temperatures was found consistently for all 
the runs owing to the thermocouple position. However, 
the interpretation of the DTA results for multicomponent 
alloys was not always clear. For example, in the 
Al-Cu-Si-Mg case four peaks were obtained, one peak 
corresponding to each of the intermetallic compounds 
present in the microstructure. The temperatures of 
these peaks did not correspond to the melting temperatures 
of the compounds that were previously known. 

From an experimental point of view the 
characterization of congruently melting alloys as heat 
storage materials poses no great problems. The commercial 
calorimeters performed well within the expected limits, 
uncertainties of 5% for DTA and 3% for DSC. Conversion 
of heat storage densities from a mass to a volume basis 
requires improved density measurements. The present 
uncertainty in the room-temperatures densities is due mainly 
to the porosity of the alloys. The use of the induction 
furnace, with its vigorous stirring yielded alloys with low 
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porosity. Measurements of densities at the transition 
temperatures are being performed in a related study (20) . 

Although the DCS gave reliable heat capacities for 
alloys that melt below 700°C, for the higher temperature 
alloys the measurement of heat capacities of the liquid 
should be done by some method other than DTA. Vaporization 
losses and reaction of the vapors with the Pt cups are 
likely to be responsible for the failure of DTA to yield 
acceptable heat capacities. The vaporization problem 
does not affect the measurements of heat of fusion 
appreciably because this shift occurs to the baseline as 
a whole. However, there may be a residual effect that 
cannot be determined from the present results. 

As reported previously (19) , the heats of trans- 
formation are systematically well below the calculated 
values as summarized in Table 5. The predictions of the 
ideal model deviate to a large extent from the true 
eutectic temperatures, but the compositions are a useful 
approximation. The regular solution model improves the 
values obtained for the temperatures, but not appreciably 
those for the compositions. In the ternary case there is 
no improvement whatsoever, due to the assumptions that 
had to be made in estimating the ternary activities 
from binary data. 


Further alloying of the promising binary 
eutectics lowers the transformation temperature as 
expected, but it does not necessarily improve the heat 
storage density as expected. The deviation of the 
actual values from those calculated is erratic and 
unpredictable. The main reason for this variation 
probably is the different degrees of ordering that take 
place in the liquids. Since the presence of stable 
intermetallic compounds was sought to increase the heat of 
transformation, the probability of ordering in the liquid 
is high. To take this complication into account 
additional thermodynamic data are needed. 

The fact that in many cases the addition of a third 
component did not improve the heat of transformation as 
expected suggests that ordering in the liquid is more 
likely for multicomponent systems than for binary systems. 
Statistically it is more likely that some pair of elements 
will exhibit strong affinity in the liquid and therefore 
contribute to short range ordering. 

The practical implication of these conclusions is 
that it will be difficult to find alloys among these metals 
that have storage densities substantially greater than the 
ones presently reported, although some aspects of their 
performance such as volume change during melting may 
continue to improve. Investigation of quaternary alloys 


may put metals into a more competitive position for 
heat storage at lower temperatures. Further alloying 
may also yield a better accommodation of the melting 
temperatures to the sources. 
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APPENDIX C 


APPENDIX D 


VOLUME CHANGE DURING MELTING 
OF EUTECTIC ALLOYS 


BRUCE PREGGER 
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Radiation Source 

In previous work the source of radiation used 
was an X-ray tube powered by a Philips XRG-3000 high 
voltage generator (Appendix C) . An Ag tube was used to 
provide AgKa X-rays at 22Kev. The high intensity of this 
source necessitated the use of a graphite monochromator 
and Zr foil filters. Results for the volume expansion of 
Al, Al-0.12Si, and A1-0.18CU were obtained with this 
system that agreed with previously reported data. However, 
several characteristics of this source were undesirable. 
Variations in the observed intensity were systematic in 
nature and were attributed to a real variation in the 
output of the X-ray tube over the course of an 
experimental run. These variations were dependent on the 
change of the rate of flow of cooling water to the generator, 
and possibly also on variations in line voltage. The total 
intensity of the tube was a problem. The simple addition 
of absorbing Zr foils was not satisfactory, because high 
energy continuum radiation was relatively unaffected and 
remained to saturate the detector. A monochromator 
eliminated the unwanted radiation, though problems in beam 
alignment resulted. After alignment, however, the 
monochromator-Zr foil combination produced an intensity 
which yielded good counting statistics with an intrinsic Ge 
detector. Finally, the nature of the experiment required 
that the sample furnace be placed in close contact with 


i 
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the X-ray generator, a situation which led to the premature 
termination of experimental runs as the furnace became hot. 

A theoretical limit of 1000°C was obtainable with the 
furnace, but at about 800°C the metal base reached what 
was determined to be a maximum safe temperature. 

The use of a radioactive isotope as an X-ray source 
in place of the X-ray tube and generator offers a number 
of advantages. The output intensity of an isotope is 
constant over experimental times within the limits of 
radioactive decay statistics. Such a source is small and 
easily managed; alignment can be facilitated. The source 
can be physically separated from the furnace, allowing 
maximum temperatures. The radiation from an isotope is 
characterized by sharp energy peaks, so that the use of an 
energy dispersive counting system with such a source 
makes possible measurements without a monochromator. The 
judicious choice of an isotope, most of which have more 
than one useful energy peak, allows measurements to be made 
on materials of very different densities. From these 
considerations it was concluded that the volume change studies 
should continue with an isotope replacing the X-ray tube/ 
generator. 

The choice of an isotope was dependent upon the 
following factors--the energy peak (or peaks) had to 
optimize the volume change measurements for the greatest 


number of alloys, and the isotope had to be commercially 
available at a reasonable price. 

To optimize volume change measurements the 
equation 

(I £ /I s ) _ e " (y/p)t(p s“ p Ji ) 

20 2(/N +/BKG) 

N = number of counts - background (BKG) 

2o = uncertainty/measurement 

must be maximized. This depends upon the thickness of the 
sample, a controllable parameter, the number of counts per 
measurement and background, which are somewhat adjustable, 
and on the volume change and mass absorption coefficient, 
determined by the material. 

The alloys to be measured fell into two groups 
separated by their absorption characteristics — those 
containing Cu or Zn, and those that did not. Cu and Zn 
absorb radiation much more strongly than Al, Si, Mg, or 
P, and thus an alloy with either one must be measured with 
high energy X-rays. Radiation of this energy does not 
show enough variation with density to be of use with the 
lighter alloys. Another consideration is the detector, 
as a Nal crystal and a Ge detector behave considerably 
differently at different energies. 


Calculations were done to compare the relative value 
of two promising isotopes, Cd-109 and Gd-153. Cd-109 

has peaks at 22, 25 and 88 KeV, Gd-153 at 41 and 47 KeV. 
Calculations were made for the 22 KeV Cd and the 41 KeV 
Gd peak as they would be detected with a Ge solid state 
detector, and for the 88KeV Cd peak as it would be 
detected in a Nal crystal. A series of calculations was 
performed for each alloy at each wavelength in order to 
determine the sample thickness that would give the best 
volume change precision. Values of density were measured 
from our prepared samples. Mass absorption coefficients 

were estimated by summing the mass fractions of the 

* 

coefficients of the component elements. Elemental 
absorption coefficients are given in Table 1, estimated 
alloy coefficients in Table 2, and optimum volume change 
precisions with corresponding sample thicknesses in Table 3. 

From Table 3, it can be seen that the 41 KeV Gd peak 
with Ge detector gave the best precision for the dense 
Cu,Zn alloys. For the lighter alloys the 22 KeV Cd peak 
gave the best results. This peak energy could not be used 

* 

Volume change estimates were made by summing the mass 
fractions of reported elemental change values. The count 
rate was estimated through consideration of detector 
efficiencies and collimated radiation beam intensities 
produced by 5 mCi sources. Counting time was taken to be 
1000 seconds. 


Table 1 

Mass Absorption Coefficients for Cd-109 
Gd-153 (Cm 2 /gm) 


KeV 

Mg 

A1 

Si 

P 

Cu 

Zn In 

22 (Cd) 

2.12 

2.60 

2.7 

3.95 

24.38 

26.75 

41 (Gd) 

0.4* 

0.5* 

0.6* 

0.6* 

4.0* 

5.0* 

88 (Cd) 

0.17 

0.19 

0.2* 

0.2* 

0.67 

0.78 2.2 

Sources : 

"Handbook 

of Chemistry and Physics," 50th edition. 

C 

. R. C. Press, 1970; 

C. S. 

Barrett 

and T. 

B. Massalski 

II 

Structure of 

Metals , 

" 3rd 

edition, 

McGraw- 

-Hill, 


1966. Starred values are interpolated. 
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with the Cu and Zn alloys, but the 88KeV peak can, 
though with less precision than the Gd peak. It was 
decided to use the Cd-109 source. The better results 
for the light alloys would be useful for comparison to 
the results obtained for such materials with the X-ray 
tube used previously. 


APPARATUS 

Radiation Source 

The Cd-109 isotope used was in the form of a small 
disk, 0.79 cm in diameter and 0.51 cm in thickness. It 
was supplied by New England Nuclear (model NER-465 LE66) 
and had an activity of 5 mCi. The source is collimated 
into a narrow beam by having it rest at the closed end 
of a 15 cm stainless steel tube, held in place by a 
smaller diameter tube of the same material which slides 
into the first tube. The second tube has an inside 
diameter of 0.4 cm which acts as the collimator. The 
0.4 cm diameter corresponds to the active area of the disk 
source, ensuring the greatest possible radiation intensity 
for the collimated beam. These tubes containing the 
source are sunk into a lead brick, as shown in Figure 1. 

Sample Furnace 

The furnace is basically that used by Harrison 
(Appendix C) with some minor modifications. Five meters of 
Nichrome wire are wound about the outside of a 9 cm diameter 






alumina tube. This tube has two 1 cm diameter holes 
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drilled through it in order to let the X-ray beam pass. 

The tube sits upright in the center of the furnace on a 
refractory base. The sample holder is placed inside 
this tube. The furnace is designed to operate under 
vacuum. This vacuum is maintained in a stainless steel 
cover which fits onto an O-ring in the furnace base, and 
which has two windows built into it to allow the entrance 
and exit of the X-ray beam. These windows are made of 
beryllium, 0.01" thick, held in place by O-rings in 
stainless steel fittings soldered to the cover with silver 
solder. Thermocouple wires, electrical connections, and 
the vacuum line penetrate through the base of the furnace. 
The thermocouple is Chromel-Alumel . A vacuum of less than 
5 torr is realized with a small mechanical vacuum pump. 
Between the furnace windings and the vacuum cover are 
two stainless steel heat shields and 1 cm thickness of 
heavy Alundum refractory. The furnace is shown in Figure 2. 

Radiation Detectors 

Two types of radiation detectors were used, depending 
on the energy of the radiation being measured. For low 
energies, the Cd-109 22 and 25 KeV peaks, a General 
Electric 400 high purity germanium spectrometer was found 
to give the best results. The resolution of this device 
at 22 KeV was measured to be 2.3%. At this energy a Nal 
scintillation crystal had very poor resolution, and could 
not resolve the 22 and 25 KeV peaks and thus could not 
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measure truly monochromatic radiation, leading to 
erroneous results. However, as the energy increases the 
efficiency of the Ge counter drops rapidly. At 88 keV 
the counting rate of the Ge detector is at such a level 
that a Nal crystal is of more value as a monitoring 
device, even though the background count rate in the 
scintillation crystal is much greater. 

The apparatus is arranged on a nonconducting 
bench top in such a manner as to allow changing from 
the solid state to the scintillation counter. The isotope 
holder and furnace are both small and easily repositioned, 
being physically separate. 

Either detector was connected to a Tennelec TC 203BLR 
linear amplifier, from which the signal was run to a 
General Electric 700 multi-channel analyzer. 

The arrangement of the apparatus is given in 
Figure 3. 

Sample Holder 

The sample holder has three requirements. It must 
keep the sample thickness constant throughout the 
experimental temperature range; it must allow enough 
radiation through the sample to provide good counting 
statistics; it must be made of material that will not 
react with the samples at elevated temperature. Graphite 







was chosen as the sample containment material. Three 
graphite plates are placed together, the center plate 
cut out in the center to contain the sample. The outer 
plates are thinned to 0.2 cm in their centers to allow 
passage of the radiation. These three plates are 
sandwiched between stainless steel plates and pressed 
together tightly with eight screws. The steel plates 
have holes cut out of their centers to allow X-ray 
transmission. The configuration is shown in Figure 4. 
Graphite was chosen as containment material because of 
ease of fabrication and low X-ray absorption character- 
istics. The 0.4 cm total path length for X-ray absorption 
in the sample holder is calculated to transmit 90.9% of 
incident radiation at the 22KeV Cd-109 peak, and 95.8% 
at the 88 KeV peak. 

To facilitate alignment of the sample holder with the 
furnace windows, radiation source, and radiation detector, 
the holes in the stainless steel plates were threaded 
and fitted with 5 cm long tubes of the same material. 

These tubes defined the radiation path through the 
furnace winding and sample holder. 

Samples 

The aluminum and Al-0.12 Si samples were commercial 
materials. All other alloys were prepared by heating 




constituent elements and various commercial Mg, Si, 

P-metal compounds in an Ajax-Northrup induction furnace. 

The preparations were sealed in alumina crucibles to 
prevent vaporization of the volatile elements, and heated 
in an argon atmosphere to prevent oxidation. The induction 
heating provided good mixing of the alloys so that homo- 
geneous eutectic structures resulted. The alloys were 
machined into disks that would fit tightly between the 
two outer graphite plates in the sample holder, but smaller 
in diameter than the hole in the center plate to allow 
for transverse expansion on heating and melting. The 
sample thicknesses were determined by optimum volume change 
calculations (see Table 3) . 

The alloys themselves were those which showed 
promise as heat storage materials. Eutectic compositions 
in the various systems were found by Farkas, Appendix B. 

The heats of fusion were measured in the same work. 

Those alloys which demonstrated the greatest heats of 
fusion have possible commercial utility and were chosen 
for this study. Their compositions are given in Tables 
2 and 3. 

Experimental Procedure 

The alloy sample is placed inside the sample holder 
which is then screwed together tightly. The sample holder 
is positioned in the furnace so that the radiation passes 
from the source to the detector by way of the radiation 


path tubes on the sample holder. The radiation from a 
5 mCi Cd-109 source is too weak to permit the use of a 
fluorescent screen to check alignment; accuracy in this 
regard must be determined by counting with the detector. 
After satisfactory positioning of the sample holder in 
the furnace, the furnace cover is put in place and the 
vacuum pump is turned on. With the vacuum at 5 torr or 
less, counting may begin. The integral number of counts in 
the proper peak is recorded by the multi-channel analyzer 
for a period of 1000 seconds. A variable output trans- 
former which controls furnace heating is turned on and 
the furnace is allowed to heat so that the temperature, 
as recorded by the Chromel-Alumel thermocouple, is raised 
by about 100 °C. At this point the transformer is adjusted 
until a temperature equilibrium is achieved. The multi- 
channel analyzer is turned on to record for another 1000 
seconds. Changes is furnace temperature in this time 
period can be minimized with a skilled hand at the 
transformer; such variations are accounted for by averaging. 
The process is then continued, raising the temperature in 
100° intervals through the melting temperature and into 
the liquid range as far as the furnace allows. Measurements 
are also taken for 1000 seconds during the heating intervals 
for comparative purposes, in which temperatures are recorded 
at the beginning and end of the 1000 sec interval and 
averaged. 


RESULTS 


The values obtained for volume changes of the 
various alloys studied are given in Table 4. Figures 5 
and 6 show plots of Hn intensity vs. temperature for the 
various materials studied. The solid and liquid phase 
data points were separately characterized as straight lines 
by the method of least squares. It was felt that the 

scatter of the data precluded the use of higher order 
polynomials. These lines were extrapolated to the melting 
point of the material. The values of the solid and liquid 
line at the melting point were taken as the natural log 
of the incident intensity produced by the solid and liquid 
melting point densities, and as such were used in calculating 
those densities. The value of the solid phase data line 
was used at room temperature to represent the room 
temperature density; the measured room temperature 
intensity was not. Precision was gained in this manner. 

Numbers obtained for the solid and liquid 
expansion coefficients were not in agreement with 
previously measured values, nor were they to any great 
extent reproducible. This is a serious problem, for it 
can affect the fusion volume change calculation. The 
fusion volume changes calculated appear to agree with those 
that have been determined by Harrison and others (2-4) . 

This study is concerned primarily with new alloys on which 
no data of this nature are available for verification. 


Table 4 


Volume Changes on Melting of In, Al, and 
Various Eutectic Alloys 



(AV/V ) Measured 
s 

Tabulated 


In 

0.031 

0.025 

(1) 

Al 

0.09 

0.072 

(2) 



0.0769 

(3) 

Al-0 . 12Si 

0.041 

0.040 

(2) 


0.038 (4) 

Al-0.12Si-0.05Mg 0.066 

Cu-0.25P-0.14Zn 0.033 

(1) Smithells, "Handbook of Metals" 1955. 

(2) Harrison, Appendix C. 

(3) P. J. Wray, Met. Trans., 1974, vol . 5, p. 2602. 

(4) L. Mondolfo, "Aluminum Alloys," Butterworths , London, 1976. 
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Al-0.12Si had a change of volume on melting 
that was measured to be 4.1%, very close to the change 
of 4.0% given by Mondolfo (4) and 4.8% as measured by 
Harrison (2) with this method using an X-ray tube source. 

In, measured with the 88 KeV peak, gave a value 
of 3.1% for the volume change. Smithells (1) gives 
2.5%. 


The runs made on pure aluminum were not altogether 
satisfactory. High volume changes were calculated from them, around 
20% higher than the grouping of other experimental values 
(2-4). The reason for this remains unknown. 

Al-0.12Si-0.05Mg gave a measured volume change of 
6.6%. This value is probably too high, although there are 
no other measured values to check against. The presence 
of Mg has been found to be detrimental to the acquisition 
of good data when present in larger amounts in a graphite 
container; perhaps it is not helpful in this case. 

Cu-0.25P-0.14Zn is an easy-to-handle alloy. The 
volume change measured was 3.3%, which in concert with a high 
heat storage capacity and ease of handling (it was made 
successfully in a resistance furnace) makes this material 
look very promising for commercial utilization. 


DISCUSSION 


The variability of solid and liquid expansion 
coefficients demands resolution. The problem stems from 
the changing of the radiation intensity in other parts of 
the furnace than the sample. Due to the low level of 
intensity produced by the radiation source, as wide 
a beam as possible is used. This leads to blocking by 
opaque surfaces in the furnace. When the furnace heats, 
these surfaces may move in such a manner as to increase or 
decrease the theoretical incident intensity so that the 
equation for p mp - P ROOM T does not hold. To decrease 
the beam diameter to that usable with an X-ray tube 
source would decrease the intensity from the radioactive 
source to a level that would eliminate any chance of 
useful results. Experiments were done using an empty 
sample holder to test this theory of changing incident 
intensity, and what was found was that the furnace did 
change the intensity on heating, although not always by 
the same amount nor always even in the same direction. 

The p 0 - p is relatively unaffected by this problem coming 
from data of a single temperature, but if p is calculated 
wrong then the percentage volume change on melting will 
not be accurate. 

A solution to this problem involves the simultaneous 
measurement of two peaks of different energies. Opaque 
blocking accounting for a changing incident intensity 


may be eliminated between equations for the absorption 
at the two energies giving the following equation: 
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1 measured exp 
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where 


1 "adjusted = intens i t y that would be measured had 

not the incident: intensity been changed by 
opaque blocking from T' to T" 
y/p ^ 2 “ mass absorption coefficients for radiation 
of energies 1 and 2 

= measured intensity at T' with energy 1 (for 
example) 

I" . = I? or H . yielding corrections at both 

energies. 


This equation accounts for partial opaque blockage, but 
there is no solution for the opposite effect, partial 
free transmission. Thus any porosity in the sample 
cannot be mathematically eliminated. Samples prepared 
in the induction furnace have shown good density character- 
istics, without excessive porosity. 


Cd-109 has two peaks in the same anergy range: 
22 and 25 KeV. Test measurements were made on Al for 
corrections by the above equation. The 25 KeV peak is 


relatively weak, yielding poor counting statistics, and 
results were inconclusive. Gd-153 on the other hand has 
two strong peaks in the same range, at 41 and 47 KeV. 

An active Gd source could be a solution to the problem. 

Another problem was encountered when dealing 
with alloys rich in Mg. At elevated temperatures the 
Mg reacted with the graphite of the sample holder 
changing both sample thickness and sample composition. 

Such was the case for Si-0 . 35Cu-0 . 28Mg and Zn-0.40Cu-0.15Mg. 
A new material must be found to contain these alloys. 

It appears as though a Gd-153 source could provide 
better results than the present Cd-109. The source 
and experimental configuration now being used do not 
seem wholly satisfactory, especially for expansion 
coefficient determination. 
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ABSTRACT 

Within the field of energy, there is an increasing emphasis 
on the use of phase-change materials (PCM) for thermal energy 
storage, but the analytical and numerical methods available to 
describe their thermal behavior in cylindrical coordinates are 
relatively few in number. In addition, these methods are either 
mathematically too complex for engineering purposes, or simplified 
to the point of solving a pseudophase- change process far removed 
from the actual system itself. In this paper, an efficient, nu- 
merical method is presented for modeling one-dimensional heat con- 
duction (including sensible and latent heat) around an internally- 
convected cylinder of PCM. Finite-difference equations monitor 
the heat transfer within the solid and liquid regions, while Tny- 
lor-series expansions of temperature are made in the regions close 
to the interface. Inputs to the routine include temperature- 
dependent properties, initial temperature distributions, time- 
varying convection boundary conditions and time-varying inlet fluid 
temperatures. Outputs from the method include melt-line locations 
and temperature distributions vs. time, as well as charging rates 
and pipe surface temperatures. 


1. INTRODUCTION 

Many heat transfer processes include a change of phase of a 
material due to melting or solidification. Ice formation, abla- 
tion of missile skins and solidification of metal castings arc a 
few such applications. Within the fields of Thermal and Solar 
Knergy Storage, there is an increasing emphasis on the use of such 
phase-change materials (PCM) , for the collection, storage and trans- 
mission of the forms of energy which today are more costly to gen- 
erate and more actively conserved. The many applications of, and 
increased interest in, PCM's require accurate mathematical model- 
ing of the thermal beluivior of such materials. Predictions of 
melt-line locations and temperature histories throughout the medium 
are useful in: a) choosing the appropriate material for PCM con- 

tainment, b) optimum design of the PCM apparatus, and c) the con- 
trasting of several PCM's to choose those applicable for a specific 
use. Due to the inherent temperature non-linearity at the inter- 
face, mathematical modeling of the pha5;e-chanqe process via an 


analytic approach is quite complicated. As a result, solutions 
exist for systems with simple geometries and very restrictive 
boundary and initial conditions. The extension to systems in 
cylindrical coordinates further reduces the number of available 
solutions due to the added mathematical complexity of the 1/r 
term in the original heat conduction equation. However, this is 
the system geometry often used with PCM's and energy storage or 
transmission; such as that found with heat exchangers, heat pipes, 
tube flows, etc. 

Authors such as Meyer [2], Poots [3], Solomon [4], Boley [5], 
Krieth [6], JiJi 17], Shubin [8] and Sikarskie [14] have attempted 
such modeling, but required such restrictions as constant wall 
temperature, constant material properties or temperature distri- 
bution assumptions. The intricate mathematics or system-simplifi- 
cation of such methods preclude their use for engineering purposes. 
This is also evident in the numerical solutions alternately em- 
ployed to model the phase-change process. Talwar [9] and others 
reduced the problem to one of pseudomelting whereby the thermal 
properties at the interface were "smoothed out" to provide a con- 
tinuous relationship between the solid and liquid regions. Forster 
[10] considered the use of a quadratic temperature profile through 
the melt-line between two known points in an effort to compute the 
temperature gradients at the melt-front and hence its velocity. 
Murray and Landis [11] further refined Forsters model by institu- 
ting temperature profiles three nodes from the interface, but re- 
sulted in starting and terminating errors that propagated through- 
out the solution. Lazaridis [12] attempted to expand Murray and 
Landis' work to multidimensional systems, and resulted in develop- 
ing a routine both mathematically and computationally difficult to 
implement with various geometries. 

This paper outlines an efficient, numerical routine for solv- 
ing one-dimensional heat conduction (including sensible and latent 
heat) around an internally-convected cylinder of PCM, by noting 
the dominant mode of heat conduction to be in the radial direction. 
Finite-difference equations monitor the heat conduction in the re- 
gions away from the fusion front, while Taylor-series expansions 
of temperature are made near the melt-line, in such a manner as to 
continually trace its propagation throughout each node, while ad- 
justing all nodal temperatures. The versatility of the method is 
reflected in the fact that temperature-dependent properties and 
initial temperature distributions may be included as well as time- 
varying boundary conditions. Time-dependent inlet fluid tempera- 
tures can also be handled which represents the more realistic 
condition encountered with energy storage and transfer systems 


2 . PROBLEM STATEMENT 

The heat exchanger geometry under consideration is depicted 
in Figures 1 and 2. Figure 1 is an axial view of a cylinder sec- 
tion of length containing a PCM partitioned into "N" equal- 

thickness, concentric, annular regions or nodes, with a represen- 
tative melt-line location. Heat convection occurs at the inner 
radius Rj_, and, for the analysis, one-dimensional solidification 
occurs radially outward. The outer radius \\ 2 > is an insulated 
boundary, but may also be subjected to convection if required. As 


the convecting fluid tra- 
vels through the cylinder 
section an energy balance, 
described in Section 4.5, 
is performed on its bulk 
temperature to arrive at an 
increased outlet tempera- 
ture. This in turn becomes 
the inlet flow temperature 
for the next downstream 
section. This process is 
repeated for "M" sections 
until the desired pipe 
length of "Mi" is attained, 
as depicted in Figure 2. 

In this manner, by taking 
a sufficiently large number 1 
of axial segments, the heat 
exchanger can be effective- 
ly modeled as a two-dimen- 
sional system with the one- 
dimensional radial heat 
conduction mode dominating 
the axial conduction mode. 



Figure 1. Axial view of heat-exchanger 

with representative concentric 

3. PROBLEM SOLUTION rings (nodes) of PCM, depicting 

a particular melt-line location 

The only requirements 
of the following analysis 

are that the fusion temperature is constant, the density is con- 
stant in both phases, and that natural convection effects are neg- 
ligible. Use of the foregoing restrictions is acceptable for many 
engineering purposes, while exclusion of the same would result in 
a more complicated system modeling. As previously mentioned, only 
radial heat conduction will be considered. 


The problem analysis is facilitated by the introduction at 



SECTION I SECTION 2 


SECTION M 


Figure 2. Horizontal view of heat exchanger of length "MS," contain 
ing a PCM between an internally-convected pipe and an 
insulated outer one. 




this point, of the governing equations used throughout the phase- 
change medium: 
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where subscripts 2 d 2 denote the solid and liquid regions re- 
spectively. Eqs. la, b are the field equations in the solid and 
liquid regions. Eq. 2 is the well-known Stefan Condition result- 
ing from an energy balance at the interface, and Eq. 3 is the con- 
dition of constant fusion temperature. As an added convenience, 
the following non-dimensionalizing parameters are introduced into 
the above equations: 


£ e r/R ; r E a ± t/R 2 ; e(i) = S(t)/R ; AF, e Ar/R ; K . * 2 /k l 


y e E o 2 /a 1 ; X = L/C, ; = 6r/R ; ^ = R n ; f, M e R 
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where R is a chosen reference length. The resulting equations, in 
which all variables except temperature are non-dimensionalized, ap- 
pear as : 
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Figure 3 illustrates the cross-sectional, radial view of a 
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Figure 3. Radial Temperature distribution near the interface, with- 
in a particular pipe section. 


representative pipe section "i", 1 < i < M, undergoing solidifica- 
tion at some intermediate time "t", where the melt-line has just 
traveled to the distance ei(x). The node in the solid region 
closest to the interface is denoted as "q", while "q + 1" repre- 
sents the similar node in the liquid region. The nodal spacing, 
A£, is defined as (£n “ £i)/(N - 1) while is the time-varying 
distance between node q and ej.. At this particular time x, only 
the current location c i is known, as well as the "old" (x - Ax) 
nodal temperatures Tj^j; i = 1, M; j = 1, N. The calculation of 
"current" temperatures and future melt-line positions requires the 
successive execution of the following computational logic: 

a) Determine all nodal temperatures except Tj,,q and Tj, ,q+i, 
at time x via the finite-difference equations, 

b) Compute the current temperature qradients at the inter- 
face, 3 Tt/ 9£;| = A^ and STo/'d^ | = and, 

c i e i 

c) Determine the current values for temperatures Tj^q and 
T i,q+1' as we H as fc he melt-front velocity de^/dx? 

all of which are outlined in the next three sub-sections. 



3.1. Finite-Difference Equations 

Temperatures at all nodes i,j where 1 < i < H and 1 _< j < q-1. 



q + 2 j < N are computed from the usual finite-difference approx- 
imation to the original heat conduction equation (4) , 



where T^j denotes the new value for j after the time increment 
At, and y - 1 when j is within £he soliaregion or y = Yj, when j 
is within the liquid region. Ti f j is solved from (6) as. 
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yielding 

the first 

familiar 

stability criterion 


At < AC 2 / 2Y (8) 


Two modifications of (7) are required for the two boundary condi- 
tions specified in the problem description. 

Insulated boundary . The outer boundary temperatures are de- 
termined from (7) by setting T^j+i = Tj. # j-l» thus 




T 

i,N-l 


(9) 


Convection boundary . Using the fictitious node concept posi- 
ted by Ozi$ik [13], the nodes Ti,j-1 for j = 1 are defined as Tfi c ., 
and their values are arrived at via an energy balance at the con- 1 
vective boundary 



( 10 ) 


where Tf. is 


the bulk fluid temperature at location "i". 


Solving 


for the fictitious nodal temperatures in (10) yields. 


’*lo i " Bl<T t . - T.p + T i 2 (11) 


where the Biot number Bi = h R/k^ Inserting (11) into (7) yields 
the relation 



From (12) arises a more restrictive time stability criterion, 
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to be used in place of (8), and a second, milder, geometric criter- 
ion from the coefficient of the third term in that 


2 1 

>0 or AS < 2£ 

AS S x 


which implies that the mesh size should be sufficiently smaller 
than the cylinder inner radius. Thus, equations (7) , (9) and (12) 
are used to determine all "current" nodal temperatures except Tj^q 

and T i,q+1* 


3.2. Interface Region 


At this point, Taylor series expansions of temperature are 
imposed from the melt- line to the two known temperatures, Tj^g + 2 
and Tj^q-.]., in the form 
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generating the following two equations: 
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containing the unknowns A^, and the two second derivatives. In 
addition, (17) is valid at the melt-line for all times. 


dT k 

3T k 

0 — — 

3T k 

dx 

e . 3x 

G i + 



k = 1, 2 


(17) 


where the first quantity on the right hand side is merely the heat 


conduction 

equation 

(4) 

. Thus, (17) 

gives rise to two more equa- 

tions 
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with the introduction of another unknown de^/dx, which is the Ste- 
fan Condition (5) rewritten below as, 


de . 1 

— - = - (A. - KB. ) (19) 

dx X 1 1 


Hence, by solving (15) and (16) 
and substituting (19) into (18) 
the two unknown gradients. 


for the second derivative terms 
yields t^o coupled equations with 



( 20 ) 


( 21 ) 


where = A£ + <5Ci» d 2 = 2A£ - These two equations are then 

applied at each "ith" location, where A^ and are found via an 
iterative process, such as a Newton- Raphson routine. 


3.3. Determination of Tj^g and Tj^g+i 

With the known temperature gradients, (15) and (16) can be 
solved for the two second derivatives, in which case everything 
is known about the Taylor-series expansions made near the inter- 
face. Of further significance is the fact that such expansions 
imply a quadratic temperature distribution throughout the inter- 
face regions which is evident in the similarity between (14) and 


T(C) = a + bC + c£ 2 (22) 

where a = T m , b = A^ or Bi and c = 3 2 T 1 /3C 2 | c . or 3 2 T;>/3£; 2 | r . 

depending upon the region. Thus, temperatures at 1 nodes "q" and 1 
"q 4- 1" can be calculated from (22) by using £ = - <$Ci and £ = A£ 

- respectively. 

Lastly, the new melt-front velocity at each "i" location is 
computed from (19) . In general, the future melt-line position 
within each pipe section is determined by multiplying the melt-line 
velocity times the time-discretization parameter Ax (this repre- 
sents an inherent limitation in the finite-difference method where 
the values calculated are to remain constant over the time interval 
Ax). However, due to the cylindrical nature of the problem, i.e. 
the increase in surface area faced by a radinlly expanding inter- 
face, the Stefan Condition which heretofore was appropriately 
utilized in the mathematical derivations at the melt-front, must 
be interpreted differently when dealing with a time-increment ana- 
lysis. The energy balance at the interface, assuming unit depth, 
is 


3T. 3T 2 pLV 
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(23) 


where V is the volume of material undergoing a change of phase dur- 
ing the time increment, or "^t-rAt - r£) # ariv. r t is the radial 
melt-line location at time "t". In non-dimensional form, (23) re- 


duces to 


A i " KB i “ * 


’t+At 


which is similar to (19) multiplied by an adjustment factor that 
reduces to unity as At is made relatively small. With the newly 
computed value of melt-line location at each "ith" section, the 
process outlined in 3.1 to 3.3 is repeated anew, and in this man- 
ner the melt-line propagation and nodal temperatures are continu- 
ously monitored throughout the medium. 


3 . 4 Starting Procedure 

For the initialization analysis, the more general case of 
superheat will be taken and it is assumed that the entire medium is 
subject to a temperature distribution containing temperatures cjreat- 
er than Tm. Heat convection imposed at the inner radius success- 
ively reduces the wall temperature, at a particular- "i" location, 
to a point where change of phase is about to begin (t = 0) . At 
this time the temperature distribution throughout the liquid has 
been determined by the finite-difference equations alone. From 
this existing distribution, the first three nodal temperatures are 
used in (22) to determine a temperature distribution, from which 
3T 2 /3£ L can be determined, in a manner similar to that of [16], 


■ ( 2 T - 2 - ¥ - I T ")/ AC 


Alternately, we can solve for the second derivative term in (22) , 
and use (18b) to solve for the fusion front velocity during the 
first time increment. 


= IjLIa | _ !i\/b. 

t= 0 \ H 2 q// 1 


Although theoretically the melt-line velocity initiates at infinity, 
the above analysis provides a usable value of velocity for input to 
the finite-difference routine. This novel approach obviates the 
more difficult task of first determining the values A^, which now 
are more easily derivable from (19) as 


df i 

A. = — =• A + K B. 

1 dx X 


T = 0 


I* 


"l"th pip* section 

Figure 4. Cross-secticnal view of the heat-exchanger with a smooth- 
line approximation to the "M” melt-line locations, and 
an elemental fluid flow section, at an intermediate time 
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line locations of the "M" axial sections, as well as a representa- 
tive pipe flow section. An iterative method is employed to deter- 
mine the average bulk temperature Tf ^ , whereby successive energy 
balances are performed on the volume of fluid passing through the 
particular node "i" . During the first pass in ^he method, Tf. is 
assumed to be equal to the incoming fluid temperature Tf^ ^ Thus, 
an energy balance on the fluid element (31) yields the maximum fluid 
outlet temperature Tf^ QUt 

Q -(hA T ifl - T fi ) = mC ( ( T£ . out - T fi(in ) (31) 


Of course, the actual bulk fluid temperature Tf^ is higher than 
that of T fi # i n due to the heating of the fluid as it passes through 
the node, and is now approximated by averaging the newly computed 
values as 


Tf i,in + T *i, 


out 


| The energy balance (31) is again applied to yield a smaller value 

* of T fi f0 ut and above process is iterated until successive cal- 

culations of Tf ^ are within a desired tolerance. At this point 
| Tf i+l,in = T fi,out' and the method is carried out repeatedly until 

I 


w 



the fluid exits the pipe. The fluid temperatures are adjusted after 
each time step used in the finite-difference routine. 


3.6 Melt-Line Crossing a Node 

When a newly computed melt-line at an "ith" location indicates 
that the melt-line has crossed over node "q+1", an error arises in 
the computation of Ti,g-1 J T (formerly Ti # qU_^ T ), whereby the tem- 
perature Ti,q | x -at = 7m + is erroneously used in the finite-differ- 
ence routine to compute Ti,q_i | T (Actually, during the time-incre- 
ment, its temperature varies from T m + to T m ”, or, a value less than 
the melting temperature). To avoid this inconsistency, the time 
step is adjusted (reduced) so that the melt-line falls exactly on 
the node q+1. At this time, the Taylor-series expansions are prop- 
erly carried out, the new melt-line time location is determined by 
resuming the full time step, and the normal-procedure is followed. 


3 . 7 Terminating Procedure 

In the case of an insulated boundary at £n* the following 
routine is used as the melt-line just crosses the node An- 1- Simi- 
lar to the starting solution (3.4), the nodal temperature Tj., q +2 is 
non-existent, and the Taylor-series (14) is made only to node Aro- 
using the injulated boundary condition 



= 0 


in (14) where d = AA - 6£i» solving for the second derivation term, 
and substituting this into (18b) yields 


B 


2 

i 


±_ + ¥i 

Kd e i K / K 


0 


(32) 


which replaces (21) when the melt-line is between the nodes Aro-1 
and An* The current value of Tj^ro is then computed from (22) which 
reduces to 



dB. 


T + 
m 


(33) 


3.8 A Third Stability Criterion 

The third and most restrictive stability criterion results from 
an analysis in Appendix 1, (adapted from 116)), whereby the start- 


ing solution used In Sec. 3.4 to generate nodal temperatures is 
meshed with the finite-difference routine yielding the same nodal 
temperatures. The resulting criterion is that the thermal parame 
tors and node size should be chosen such that B AT, > 1.0 or 


h 

— Ar > 1 (34) 


Violation of (34) results in starting solutions that yield fluctu- 
ating temperatures and melt-line velocities. 

This condition of course will not arise from the use of a con- 
stant wall temperature at the inner radius. The problem solution 
with this boundary condition can be attained by the method outlined 
in Appendix II, or by the use of convection with a very high value 
for the heat-transfer coefficient, in which case (34) would again 
have to be satisfied and result in an extremely small time-step. 


4. RESULTS AND CONCLUSIONS 

The following representative program outputs were generated 
from the thermal property data of two types of phase-change mater- 
ials currently under investigation, namely, salts and metal alloys. 
The purpose of this is twofold: first, to provide usable data con- 

cerning each PCM, and second, to emphasize the ability of the pres- 
ent program in contrasting several PCM's to determine their applica- 
bility. A Burroughs Corporation B-7700 series computer was util- 
ized to illustrate the behavior of three chosen PCM’s, the thermal 
properties of which appear in Table 1. The heat-exchanger system 
is similar in function to that used in [15] where a workable config- 
uration was developed for use in a Solar Brayton Power Plant con- 
cept, and its dimensions appear in Table 2, where an effort was 
made to contrast the materials by choosing equal governing parame- 
ters (e.g. constant h and constant temperature difference between 
the melting and inlet convecting temperatures) . The assumptions used 
in producing the results are listed in Table 3, where only 1 through 
4 are required by the foregoing analysis. 

Figures 5-7 depict the melt-line locations within the heat- 
exchanger at different times, for the salt and metal-alloy PCM's. 
Actually, each interface location represents a contour surrounding 
the internally-convected pipe, so that in effect the three-dimen- 
sional shape of the melt-line resembles a cone in Figure 5, or a 
trumpet end in Figures 6 & 7. Due to the higher conductivity of 
the alloys, more heat is convected within the upstream pipe sections 
than that by the salt, hence, the contours are more curved. The 
slight curving of the fusion rronts with time in Figure 5 is attri- 
butable to the decreasing fluid temperatures within the downstream 
pipe section; as less heat is conducted into the upstream sections, 
its temperature rise is smaller. The lower thermal conductivity of 
the salt yields an expected increase in the time required for solid- 
ification. The similarity of the metal alloys is reflected in the 
similar interface shapes, with the Mg-Si alloy requiring more time 



Latent 

Heat 

(kj/kgj 

Specific 

Heat 

(kj/kg/°C) 

598.0 

.80 

515.0 

.87 

1,212.0 

.87 


Thermal 

Conductivity 

(w/m/*0 


4.3 


75.0 


75.0 


Density 
(kg/m. 1 ) 

Melting 

Temperature 

CC) 

2,339.0 

640.0 

2,250.0 

579.0 

2,000.0 

947.0 


Table 1. Thermal property data for three representative phase-change 
materials . 



Tabic 2. Heat exchanger data. 


1. Neglect natural convection 

2. Constant melting temperature 

3. Density is a constant in both phases 

4. Neglect axial conduction within PCM and working fluid 

5. Constant material properties 

6. PCM initially superheated 

7. No crystallization at interface 

8. Heat convection at inner radius 

9. Insulated outer radius 

10. Linear bulk- temperature distribution within each pipe section 

11. Neglect pipe thermal capacity 

12. Constant heat-transfer coefficient and inlet fluid temperature 


Table 3. Assumptions used to generate results in Figures 5 - 12 . 


to solidify due to its higher latent heat of fusion. Due to the 
small time step (high h value required from (34)) the melt-line 
propagation only up to about 5 cm. was plotted to facilitate compu- 
ter generation of the results; however with greater access to com- 
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line location at six different tines within the 
exchancer for the salt PCM. (3iot No. * 2537.0) 



ielt-line location at seven different times within 
he heat- exchanger for the metal alloy Al-Si. 

Biot No. ■ 145.5) 





5.50 



Figure 7. Melt-line location at seven different time within 
the heat- exchanger for the metal alloy Kg-Si. 

(Biot No. s* 145.5) 
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Figure 8. Melt-line location within the heat-exchanger for times up until 
the salt PCM is almost totally solidified. The heat-transfer 
coefficient and mass flowrate were reduced to 731 w/m2- 0 c and 
1.4 gr/sec, respectively, in order to exhibit the melt line 
propagation along the insulated outer boundary. (Biot No. * 34 
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puter process time, the entire solidification region is similarly 
produced. As an example of this, Figure 8 illustrates almost the 
entire region of the salt PCM being solidified, via the reduction 
of both the heat- transfer coefficient and mass- flowrate. The inter- 
face movement along the insulated boundary is shown, as well as the 
increase in solidification time. Also of interest, is the effect 
of mass flowrate on the interface contours. Figures S and 8 indi- 
cate that the mass flowrate can be reduced to the point where even 
the salt, with its low thermal conductivity, can transfer much of 
its energy within the upstream nodes, in a manner similar to that 
of the metal alloys. 

The efficacy of the alloys over the salt in charging a heat 
storage system is shown in Figure 9 where the area under the curves 
represents the total energy transferred after a certain time, and 
is a .measure of the effectiveness of the heat exchanger. The values 
for Q are determined by the mass flow rate, m, specific heat of the 
fluid, Cf, and the temperature difference between the inlet and 
outlet fluid temperatures. Of course Q max occurs when the tempera- 
ture difference is the greatest (i.e. initial temperature of the 
medium minus the inlet fluid temperature) . 

Figures 10 and 11 exhibit the radial temperatures within the 
first pipe section at various times, for the salt and Al-Si data 
respectively. The dissipation of heat with time throughout the 
superheated region is evident, as well as the temperature non-lin- 
earity at the interface. The convecting surface temperature is also 
depicted as decreasing with time. The fact that the inside surface 
temperature of the pipe is lower for the salt than the alloy, and 
not vice versa as one would be lead to believe, is due to the heat 
convection value chosen for the analysis. Similar to the resistance- 
network concept, the convection resistance at the wall is greater 
than that encountered within the alloy medium, but of the same order 
of magnitude as that within the salt. Hence, the rate-determining 
step for heat transfer out of the metal alloy, is the convection 
resistance (unlike the salt) , and if this were made smaller (by in- 
creasing h) then the surface temperature would decrease. In a 
sense, the heat-transfer efficiency of the alloys is even greater 
than that shown. Figure 12 provides an axial composite of radial 
temperature distributions throughout the partially solidified metal 
alloy Al-Si at a time of 37 seconds. The inlet and outlet pipe 
sections are 1 and 11, and every other section is included. The up- 
stream nodes are undergoing solidification, while the downstream 
sections are approaching the melting temperature. 
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NOMENCLATURE- (non-dimensionalized parameters appear in parenthesis) 

A^, temperature gradient at interface in solid region 
temperature gradient at interface in liquid region 
Bi Biot number , h R/ki 

Ci specific heat of solid phase 

C f specific heat of working fluid 

h heat-transfer coefficient 

ki thermal conductivity, solid phase 

k 2 thermal conductivity, liquid region 

K ratio of kl/k2 

L latent heat of PCM 

Z axial length of pipe section 

m mass flowrate 

M number of pipe sections 

N number of radial nodes 

q node in solid region closest to interface 
q+1 node in liquid region closest to interface 

R. reference length 

Rl ( C i ) pipe inner radius 

R 2 (£ N ) pipe outer radius 

r(£) radial coordinate 

Rj(5j) radial location of node 

6r^(S£^) distance between node q and melt-line 

Ar(A£) nodal spacing, Rj-Rj/N-l 

S. (e^) radial distance to interface 
t(x) time 

At (A t) time increment 

T. . fictitious node 

flCj[ 

T, fluid temperature 

r i 

T fusion temperature 

m 

T w wall temperature 

T. . nodal temperature; l<i<M;l<j<N 
1 • J — — ““ 

thermal diffusivity, solid region 
c »2 thermal diffusivity, liquid region 

p density of PCM 

A non-dimensionalized parameter, L/C^ 

Y equals within liquid region; equals 1 within solid region 
Y^ ratio of thermal dif fusivities a 2 /a^ 


APPENDIX I. DETERMINATION OF STABILITY CRITERION 

An energy-balance at the convection boundary during the onset 
of solidification yields, 


3T 

H 


Bi(T x - T f ) 


(1) 


where the subscript "i" has been dropped for simplicity, and (25) 
provides an expression for 3T/3C • Solving for the temper- 
ature at the boundary yields ^1, t«0 


4T 


2 - T 3 + 2BiAC T f 
3 + 2Bi AC 


( 2 ) 


Now after an infinitesimal time step At, (2) should also hold for 
the new values of T 3 , T 2 and T 3 , denoted by primes as. 


4T 


2 - T 3 + 2BiAC T f 
3 + 2BiAC 


(3) 


but T 2 and T 3 are found from the finite-difference equation (7) , as 
T 2 » T 2 a + b(T 3 + T^) + c(T 3 - T 3 > (4a) 


T 3 - T 3 a + b(T 4 + T 2 ) + c(T 4 - T 2 ) (4b) 

where a = 1 - 2yAt/AC 2 , b ■ yAt/AC 2 , c = Aty/2CAC 
Substituting (4a, b) into (3) yields 

T^ *» 1/d [T^ (4b - 4c) + T 2 (4a - b + c) - T 3 (a - 4b - 4c) 

+ T 4 (-b - c) + 2ACBiT f ) (5) 

where d = 2A£Bi + 3, and 

(5) can be written in simpler form as 

T^ - + 3T 2 - nT 3 - 4>T 4 + g>) (6) 


where 



At , / 1 

1\ 

0 * 

4 ( 

— ) 


AC \AC 

2C/ 


9yAT 

Aty 

6 - 

4 + 



AC 2 

2CAC 


6yAt 

4Aty 

n - 

1 



AC 2 

2CAC 


YAt Aty 


♦ ■ 

+ 



AC 2 2CAC 


(1) B 

2ACBiT f 



(7a) 

(7b) 

(7c) 

(7d) 

(7e) 


The error obtained in (6) is a result of individual errors in all 
the computed temperatures, hence 


At| = (aATj^ + BAT 2 - nAT 3 - (j»AT 4 )/d 


( 8 ) 


and the maximum error obtainable from (8) occurs when all the tem- 
perature errors have the same sign, hence. 


at' = 6, = I | o | j AT, | + |B ! i AT 2 I + In | |AT 3 | + U||AT 4 |l/d 
• L max x x * 

(9) 


Since the temperature T, changes the most due to the heat convec- 
tion at the boundary, its error should be the largest of all asso- 
ciated errors, hence, 


at! = «' < «,/d + 1 8 1 + In! + 1 * 1 ) 

1 max A 


(10) 


If 6! in (10) is less than ( the right hand side of the equation, then 
the successive values of T^ at later times is assured to contain 
an error converging to zero, and in a sense a bound is determined 


for the maximum tolerable error used in the finite-difference equa- 
tions. Fran (7d) $ is always positive, and (7a) contains the second 
stability criterion where AC << 2Ci so a is always positive. Thus 
(10) alters to 


6 X 1 6 i/ d + |$ | + | n ( 3 (11) 


and stability is achieved when 5^ < 6 ^ which implies 


o + ♦ + |6| + |n| < d (12) 


Rewriting (7) with v « Ary /AC 2 yields 


AC 

o = 4v [1 ) (13a) 

2C 


vAC 

e * 4 - 9v + (13b) 

2C 


2ACv 

n = 1 - 6v - ll3c) 

C 


vAC 

♦ ■ v + (13d) 

2C 


Case 1 . n is positive, which implies v(6 + 2AC/C) < 1 or v < 
1/6 in which case then 6 is surely positive and (12) is then sim- 
plified to 


1 - 5v < ACBi or ACBi > 1 (14) 


which would satisfy the above criteria. 

Case 2 . 6 is positive, which implies v < 4/9, hence the abso- 
lute value of n would be 6v + 2ACv/C - 1. After plugging into (12) 
and simplifying, this leads to v £ ACBi which is satisfied also 
by BiAC >1- 


APPENDIX II. CONSTANT WALL TEMPERATURE ANALYSIS 


The constant-wall temperature boundary condition represents 
a more complicated and paradoxical condition for use in phase- 
change problem initiation where the medium is initially at a tem- 
perature different from the fusion temperature (i.e. superheated). 
Essentially, this condition states that at time equal to zero, the 
wall temperature is at a temperature below the fusion point. This 
ignores that fact that the wall temperature is actually at the 
initial temperature at x = 0, and also that a finite amount of time 
is required to reduce it to the melting temperature, while also 
allowing a non-constant temperature distribution to establish with- 
in the superheated region. Thus, any attempt to model such a pro- 
cess is at best, an approximation to an unrealistic condition. 

One approach is to choose an initial time step significantly 
smaller than the one computed from Equation 8, Section 3.1, and 
perform an energy-balance on the small volume of material undergo- 
ing phase-change during that time increment. 


/ klA(T„ - V j 


At = pA6r 


( L + c l' T initial - V) 
and the non-dimensional distance is solved from (1) as 


( 1 ) 


6r 


6 ? = 


^ / ^1 ( T m ~ T w> \ 

* ypL + PC, <T initlal - T m ) ) 


1/2 


where the cylindrical volume effects are considered negligible for 
such a thin ring of material. Following the establishment of an 
initial 6£, the remaining procedure is similar to that of sections 
3.1 to 3.3, with the simplification that Ti in Equation (15), 

Section 3 .2 for i = 2 (or T^ _ when the interface is between the 
first two nodes) is known as the constant T w . Lastly, had the ini- 
tial temperature distribution been something other than a constant, 
such as a quadratic or parabolic, the starting solution outlined 
in 3.4 would apply just as well, and obviate the inherent diffi- 
culty in determining the solid temperature gradients, A^ , first. 
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I 

INTRODUCTION 

l Materials which undergo a change of phase are being intensely 

| researched for applications in the fields of Thermal Energy 

j Storage and Transfer. The high energy-density and desirable 

| thermal characteristics of these phase-change materials (PCM's) 

make them attractive candidates for use in such areas as 
solar energy collection and peak- load leveling [20] . Along with 
\ this increased interest in PCM's is the need for accurate 

modeling of the thermal behavior of such materials. Due to the 
i inherent temperature nonlinearity at the melt-line, the analyt- 

[ ical solutions available for solving the phase-change process 

> are relatively few in number, and require restrictive boundary 

or initial conditions [1-8] . The numerical approaches 
alternately employed for the problem solution, although some- 
what more applicable, also require restrictive conditions which 
I many times reduces the analysis to that of a pseudophase-change 

| process far removed from the actual system itself. Extension 

| to two-dimensional cyclindrical coordinates further reduces the 

number of available solutions in the literature [9-15] . Talwar 
and Dilpare [9] appears to be the first attempt at modeling 
^ such a phase-change problem. 

The method presented in this paper evolved from that initially 
posited by Forster [10] , refined by Murray [11] , and subse- 
quently utilized by Lazaridis [12] , resulting in a complicated 
routine difficult to implement in cylindrical coordinates. The 
feasibility of this solution was shown by Farouk [19] , but the 
required clarification and inherent simplicity was demonstrated 
by Harding [16] and Labdon [17] . 

The proposed solution method represents an efficient numerical 
routine for modeling two-dimensional heat conduction (axial 
and radial) within an annulus of phase-change material. 
Finite-difference equations monitor the heat conduction in the 
* regions away from the interface, while truncated Taylor-series 

expansions of temperature are made near the interface, so as 
to couple the thermal behavior between the solid and liquid 
regions. In this manner, the exact melt-line locations and 
temperature histories are continually determined as the 



interface propagates throughout the medium. Inputs include time 
and spatial?* dependent boundary conditions, temperature-dependent 
properties and initial temperature distributions, while outputs 
include interface positions and temperature distributions as 
functions of time. 

PROBLEM STATEMENT 

The cylindrical geometry / \ 

under consideration is / ^ \ 

depicted in Figures 1 / \ 

and 2. Figure 1 is an / \ 

axial view of a hollow j / /jy \ 

cylinder of length / / /fry \ \Va\ \ 

containing a PCM parti- / / /// /\. / \ YA \ \ 

tioned into 'N‘ equal- / If/// ! i Ylv i In 1 1 

thickness, concentric, I I'll 00* * * I |'| * * * I 

annular regions or nodes, 1 l V /\ I PI I I 
showing a particular \ \ \ / \* 2 J 111 J I 

melt- line location. For y \ Y\\ W, \ y /// J J 
the following analysis, \ \ \X V /// / j 

heat convection is \ / / 

imposed on the inner \ \ / 

cylinder Rj_, while two- \ / 

dimensional solidifica- N. \ y 

tion occurs radially X. \ 

outward toward the ^' s v ls ^ 

insulated boundary R 2 - 

Figure 1. Axial view of the cylinder 
with representative concentric rings 
(nodes) of PCM, depicting a particular 
melt- line location. 


The concatenation of 'M' cylinder sections results in a heat 
exchanger of desired length 'd', as illustrated in Figure 2. An 
energy balance is performed upon the bulk-temperature of the 



SECTION 


SECTION 2 


SECTION M 


Figure 2. Horizontal view of the hollow cylinder of length 'd 
where d ■ M£, containing a PCM between an internally convected 
pipe and an insulated outer one. 



convecting fluid as it travels through each cylinder section. In 
this manner , the melt-line contour is continually traced through- 
out the phase-change medium, as illustrated in Figure 3, which 
also shows the two-dimensional mesh network. 

PROBLEM SOLUTION 

The basic assumptions for the analysis are that the fusion temper- 
ature is constant, no crystallization occurs at the interface, 
the density is constant in both phases and that natural convection 
effects are neglected. The governing equations and conditions used 
throughout the medium are: 


3 i 2 T. 3 2 T, 

k k 

3r 2 3z 2 


1 3T, 

+ i 

r 3r 


1 3T, 


t > 0 


(l,ab) 


f k Kl 


r=S(t) , t > 0 


i th ptp« MCtion 

Figure 3. Cross-sectional view of the cylii Vjr with a smooth- 
line approximation to the 'M+l', radial melt-line locations, and 
an elemental fluid flow section, at an intermediate time 't'. 
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T ■ 

)l 

T £ (z,r) t ■ 0 
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where k = s or l depending upon the solid or liquid region. The 
following non-dimensionalizing parameters are introduced into 
the above equation. 


£ = r/R ; t = a s t/R 2 ; e(x) 5 S(t)/R ; AS = A r/R ; K 5 k t A s 
Y^ = cx^/ a g ; X = VC g ; 6£ = 6r/R ; ^ ^ ^ ^ E Z//R 

(where R is a chosen reference length) generating the following 
equation in which all variables but temperature are non- 

dimensional ized : 



3 2 T r 1 

3T k 

1 

3T k 

i * i k = s 


+ — 

= 

— 

1 ■ 


H’ 

3£ 2 S 

3£ 

Y 

3x 

e 1 5 i C N » k = i 


1 + 





S “ e (x) 


(7a, b) 


( 8 ) 


3T 

— = Bi(T -T f ) 
3 r w r 


S - c x , T > 0 (9) 




. 3T k 


= d \ 

3£ 


a; 

ii 

o 

3C 

H 

3 

S - e(t) 

» T i " 1 


( 10 ) 


( 11 , 12 ) 


'*/Kr 


5 


Method 1 - One- Dimensional Taylor Series Expansion 

Figure 4 describes the temperature distribution within the inter- 
face region for a particular cylinder section. The axial temper- 
ature distributions in the ± c direction (not shown) account 
for the two-dimensional movement of the interface. At this point, 
all that is known are the 'current' melt-line locations at time t, 
and the 'old' nodal temperatures at time t-At. The object of the 
method is to determine the current nodal temperatures and new 
melt-line velocity, which requires the successful execution of the 
following computational logic : 


Step 1 - Determine Melt-Line Slopes . The more arduous task 
of monitoring fusion- front propagation along two orthogonal 
directions is superseded by the method of determining the radial 
projection of the normal velocity gradient, shown graphically in 
Figure 5. The mathematical derivation of this procedure appears 

in (13], the result of which is a factor [l + (|^-) J = B 2 to be 

multiplied by the normal gradient, that allows continuous moni- 
toring of the melt-line gradient in the radial direction. The 


term represents the slope of the melt-line in the c direction, 
and is determined at each 'i' location as (e^ +1 -e^_ 1 )/(2Ai;) • 


T 



Figure 4. Radial temperature distribution near the interface, 
within a particular pipe section. 
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Figure 5. Particular melt-line location depicting radial projec- 
tion of normal directional gradient , and distances used for 
localized melt-line slope determination. 


Step 2 - Nodal Temperatures . The current nodal temperatures 
' j where 1 _< i < M ; 1 < j < q - 1 and cr + 2 <_ j <_ N are computed 
via J the finite-diTference equation 


T ! . - t. (l - _ 2lAl\ + m / T . ♦ T. . \ 

1,3 1,3 \ AC 2 A ; 2 / A ; 2 \ 1+1,3 1 ' 1 ' 3 / 

( T i.:+1 - T i. J-l) + ( T i. j+ 1 + T i, J-l) (l3 > 


Y At 

2U5 


with an associated time stability criterion, 


At 


< i /-L- ♦ -L- - ♦ bi \ _1 

2y lAC 2 At 2 25 A5 1 


(14) 


determined from the convection boundary [18] . 
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Step 3 - Determine Axial Temperatures and Gradients . The 
four axial, inter-node temperature* shown in Figure 6 (dashed 
circles) are required (at each location) for determining 
temperature gradients and second-derivatives at the melt-line. 
Each temperature is interpolated from a three-point distribu- 
tion determined from the radial temperature profiles originating 
from the interface (solid circles) . Using these inter-node 
temperatures and the melting temperature, two quadratic profiles 
in the form 


T (c) *a+b?+cc a 


(15) 


are utilized, both of which originate at the interface and extend 
in the i c directions (horizontal dashed ellipsis of Figure 6) . 
The two second derivatives of temperature in the axial direction 
and then determined at the melt-line from (15) and are referred 
to as Zjj, ZJ. 


Step 4 - Determine Radial Temperature Gradients . Radial 

Taylor-series expansions of temperature are then applied from 

each melt- line location in the ± £ directions (vertical 

dashed ellipsis of Figure 5) to the known value of nodal 

temperature T. „ . and T. in the form: 

i»q-l i,q+2 



U 




(£ - H ) 2 3 2 T k 


2 



(16a, b) 


£ 


Figure 6. Local tem- 
peratures and distri- 
butions required for 
temperature derivative 
calculations at a 
particular melt-line 
location. 
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In addition, (17) applies at the interface 


dT. 


dx 


3T, 


3x 


!!> 

35 


dx 


(17a, b) 


The seven equations (7a, b) , (8), (16a, b) , (17a, b) are then 
combined in a manner similar to that of [18J to solve for the 
seven unknowns, namely , the two second derivatives of temperature, 
the two radial temperature gradients, the two time derivatives 
of temperature, and the melt-line velocity, all of which are 
applied at the current interface location. The resulting two 

3T_, 3T.| 


equations are coupled in the unknowns 


and B^» 


H 


and — 


defined as 


e i 


35 


e i 


( X 2a\ 2A i v „ 

“ + H “ A i B i Ke * + ( T i,q-1 - *»)+ z s “ 0 

( 2yX Y X \ A.B.e 3 2 yA . v 

I - + — — (T - T. ) + Z" » 0 

Kd 2 K c i J K d 2 2 K \ m i,q 2 / 1 


(18) 


(19) 


Step 5 - Calculate Interface Region Nodal Temperatures . The 
similarity between (15) and (16) indicates that the Taylor-series 
expansion of three terms is identical to instituting a quadratic: 
temperature profile, and thus the foregoing analysis of the inter- 
face regions assures uniformity throughout the method. At this 
point, (16) can be used to determine the remaining nodal termper- 

atures T. „ and T. , . 

i,q i,q+l 


Step 6 - Determine New Melt-Front Velocity . Via an analysis 
in [18] the Stefan condition in cylindrical coordinates (8), which 
heretofore was appropriately utilized in the mathematical deriva- 
tions, is interpreted differently for use in a numerical routine. 
The increase in surface area faced by a radially expanding inter- 
face alters (8) by the factor below 


A, - KB, 


( :• ’) ( 


x+Ax 


+ e 


2 € . 


) 


( 20 ) 
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which reduces to unity as the time-discretization parameter 
At is made smaller. 


Method 2 - Two-Dimensional Taylor-Series Expansion 

As shown in Figure 7, 2 two-dimensional Taylor-series expansions 
of temperature are applied at the interface in the ± £ and ± c 
directions in the form: 


3T 

3T 

3 a T 

+ — 

(C-i) + — 


a; 

«. 95 

_ 1 ac 2 


(C“i) 


3 2 T 


3£ 5 


c. 

i 


(S-e^ 3 a T 

2 3 £ 3 ; 


(£-e ± ) (;-i) 


- + a(c-i) + btf-e^ + c(;-i) 2 + d(£- £i ) 2 + e(;-i) (21) 


where the constants represent the appropriate gracient or second 
derivative term in (21) , and the terms in parentheses represent 


£ 



Figure 7. Temperature field surfaces generated from two-dimensional 
Taylor-series expansions about the interface. 
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the axial and radial distance to a particular node of known 
current temperature. Thus, the selection of five known temper- 
atures in the neighborhood of the interface, in both the liquid 
and solid regions, is sufficient to determine the sought after 
temperature gradients A, cmd B, in (21). Intuitively, the 
closer the chosen nodal points to the interface, the more accuracy 
that would result in generating the two-dimensional surfaces of 
temperature, and an intricate computer program was utilized to 
determine the nodes of closest proximity to the interface. Thus, 
steps 3, 4 and 5 of Method 1 are superseded by the above proce- 
dure, which results in decreased computer process times for 
problem solutions. 

Information as to starting and terminating procedures, the 
fluid energy balance, the melt-line crossing a nodal point, and 
the use of a constant wall temperature boundary condition can be 
found in [18]. In particular, the starting procedure requires 
the following restrictive stability criterion, 


iL Ar > l.o 

k_ 


RESULTS 


The following system parameters were chosen to describe the 
various outputs from the method: 


Cylinder inner radius (cm) 

Cylinder outer radius (cm) 

Cylinder length (cm) 

Working fluid 

Working fluid specific heat (kj/kg- # C) 
Mass flow rate (gr/sec) 

Heat transfer coefficient (w/m 1 2 3 4 5 6 7 -*C) 
Initial superheat (°C) 

Working fluid temperature (°C) 


- 2.5 

- 8.0 
- 20.0 

- Helium 

- 5.2 


- 15.0 

- 5454.5 

- T + 10 

- T®* - 140 


ro 


In addition, the following assumptions were used to facilitate 
data generation only, and are not required by the foregoing 
analysis: 


1. Constant material properties 

2. PCM initially superheated 

3. Heat convection at inner radius 

4. Insulated outer radius 

5. Linear bulk- temperature distribution within each pipe 

section 

6. Neglect pipe thermal capacity 

7. Constant heat-transfer coefficient and inlet fluid 

temperature 
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Figures 8 and 9 depict melt-line contours at various times for 
the salt PCM 7CaF2-54KF-39NaF and the metal alloy PCM Al-Si , 
respectively. The former, having a lower thermal conductivity, 
yields the result of relatively straight contours and longer 
times required for solidification. Figure 9 illustrates the 
more curved interface shapes due to the greater heat- transfer 
within the upstream nodes. 

Figure 10 contains radial temperatures within the first cylinder 
section at times corresponding to those in Figure 8 for the salt 
PCM. The decreasing inner surface temperature and initial 
superheat with time are evident. 

Figure 11 represents an axial composite of several radial 
temperature distributions within the Al-Si PCM at a time of 
20 seconds following the onset of solidification. The upstream 
sections are undergoing phase-change while the initial super- 
heat is still being dissipated within the downstream ones. 
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Figure 8 - Melt-line contours at six different times within the annulus for 
the salt phase-change material 7CaF~-54KF-39NaF. 
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NOMENCLATURE (non-diraensionalized parameters appear in parentheses) 


A i 


B i 

Bi 


h 

k 


K 

L 

( 

m 

M 

N 

q 

q+i 

R 


K 1 (; i) 
r 2 (C n ) 

r(U 

w 

6 r i< 6 V 

Ar(AO 


Si(Ci) 


temperature gradient at interface In solid region 
temperature gradient at interface in liquid region 
Biot number, h R/k 

s 

specific heat of solid phase 

specific heat of working fluid 

heat-transfer coefficient 

thermal conductivity, solid phase 

thermal conductivity, liquid region 

ratio of k 0 /k 
*. s 

latent heat of PCM 

axial length of pipe section 

mass flowrate 

number of pipe sections 

number of radial nodes 

node in solid region closest to interface 

node in liquid region closest to interface 

reference length 

pipe inner radius 

pipe outer radius 

radial coordinate 

radial location of node 

distance between node q and melt-line 

nodal spacing (Rj-Rj^) / (N-l) 

radial distance to interface 
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t (t) 
At(Ax) 



T 

w 

T. ■ 

1*3 

z U) 


Z 


k 


time 

time increment 
fluid temperature 
fusion temperature 
wall temperature 

nodal temperature; l£i£M; 1 3 1 N 

axial coordinate 

second derivatives of temperature in the axial 
direction evaluated at the interface 


d^ distance from node q to melt-line 

d 2 distance from melt-line to node q + 1 

thermal diffusivity, solid region 

a, thermal diffusivity, liquid region 

2 

8^ = 1 +(||-) which is projection of normal direction onto 
radial coordinate 

p density of PCM 

X non-dimensional ized parameter, L/C^ 

Y equals Yn within liquid region; equals 1 within 

solid region 

Y^ ratio of thermal dif fusivities 
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CONTAINMENT OF HEAT STORAGE ALLOYS 


Nikhil S. Nagasw ;mi 


ABSTRACT 


The present work is directed towards finding 
optimum containment materials for several recently 
discovered heat storage alloys. These alloys include 
binary, ternary and quarternary systems of the elements 
aluminum, copper, magnesium, zinc, phosphorous and 
silicon, all fabricated or synthesized at the University 
of Delaware by Alan Riechman and Diana Farkas. The 
commercial viability of these systems hinges on their 
proving more cost effective as compared to several molten 
salt systems in their respective temperature ranges. 

Since it is anticipated that the containment system could 
prove to be as or more expensive than the storage media, 
this area of research warrants close study. 

The containment materials must satisfy several 
requirements, amongst them: good thermal stability in 

the temperature range, thermal shock resistance, high 
strength at elevated temperatures, and remain inert to 
attack by the alloy, that is, they must exhibit limited 
solubility, resistance to embrittlement, little or no 
decomposition, resistance to infiltration or liquid metal 
penetration (unless the infiltration enhances strength) and 
form a very slowly growing reaction layer that is protective. 


Experimental studies completed and in progress 
include: 

A study of wetting, reaction layer growth, and alloy 
penetration and diffusion into the containment material, 
upon exposure, without pressure at temperatures 100 to 
150 # C above the melting points of the alloys. 

Similar studies are also being done for the 
individual elements comprising these alloys. 

The experimental techniques employed in the 
present work include optical metallography, scanning 
electron microscopy and electron probe microanalysis. 

Six alloys systems are presently being studied. 

These include three alloys which show good heat storage 
properties on a volume basis and three alloys that are good 
on a weight basis. These alloys show promise for application 
in three temperature ranges: 550-625°, 700-800°, and 
950-1050°C. 

The containment materials being studied include 
two stainless steels, S.S.304 and S.S.316L, and a cheap 
silicon carbide compact, "Carbofrax Sic," from the 
Carborundum Company. 

In all the systems studied thus far, the SiC shows 
great promise as a containment material, while the stainless 
steels are attacked to a greater degree by the molten alloys. 
A conclusive report on all studies conducted will be 


available soon. 


INTRODUCTION 


An extensive literature survey has been conducted 
and various aspects of the containment problem studied. 

These include: liquid metal and solid-vapor phase 
embrittlement of metals and refractory structural 
materials, including mechanisms of embrittlement, 
specificity of embrittlement, rates of embrittlement, 
and nature of embrittlement; reaction of liquid metal- 
structural material couples and modes of decomposition, 
attack and material diffusion, surface phenomena and morpho- 
logies and effect on liquid metal interactions, infiltration 
of metals and refractory materials by liquid metals, silicon 
carbide — its preparation, properties and reactions, 
stainless steels — their properties and the effect of liquid 
metal environments, analytical procedures and theory 
behind SEM, TEM, STEM and diffraction studies and applica- 
tions to the present problem, heat storage alloy applications 
and future trends for research. 

The conclusions of the literature search are 
presented in the introduction to the thesis based on the 
present work. 


EXPERIMENTAL DETAILS 


Alloy Preparation 

The first method utilized in preparing the alloys 
was melting in a resistance furnace under argon atmospheres. 
The alloying elements were mixed in a graphite crucible 
which was closed with a tight fitting lid, incorporating 
a stirrer to facilitate mixing of the molten alloy. The 
tight fitting lid was used to minimize vaporization losses, 
particularly of the elements Zn, Mg and P. The set up is 
illustrated in Figure 1. 

The temperature of melting depended on the volatile 
elements Zn, Mg. P present and on the highest melting 
constituent present. The graphite crucible containing 
the elements was introduced into the furnace which is 
preheated to the melting temperature desired and flushed 
with argon prior to and during the melting process. The 
alloy is cooled in the argon atmosphere by switching off 
the furnace and maintaining the argon flow. 

This method proved largely unsuccessful owing to 
several reasons: 

1) When the temperature was maintained at the 
melting temperature of the highest melting constituent, 
excessive vaporization of the volatile elements occurs. 






2) When the melting temperature is held intermediate 


between the melting temperatures of the volatile and high 
melting elements , incomplete melting occurs and the alloy 
is extremely heterogeneous. 

3) Owing to the air trapped in the crucible 
within the alloying element mix, considerable oxidation of 
the alloying elements occurs, resulting in incomplete 
melting and mixing and in a porous sintered mass, even 
at temperatures higher than the melting temperature of 
the highest melting element. 

In seeking a solution to the above problem the 
following set up was utilized (Figure 2) . 

The alloying elements (in the form of chunks — 
not powders) were weighed and mixed in a graphite crucible, 
which is placed within a quartz tube as shown. The quartz 
tube is fitted with an air tight seal and evacuated at 
room temperature to remove almost all the air in the 
alloying element mix. The vacuum tight valve is then closed 
and the tube flushed with argon. In this state, the quartz 
tube placed within a set of water cooled Cu induction 
coils is heated. An argon pressure in the tube is maintained 
by regulating inlet-outlet pressures, while the graphite 
crucible and "chunky” elements provide for good coupling 
with the induction coils. Temperatures are measured by 
means of a thermocouple and an optical pyrometer. 

Following the melting, the alloy is cooled under argon. 




This system provides for quick, convenient melting 
of the alloys and minimizes vaporization losses. In 
addition the induction furnace set up provides good mixing 
of the alloying elements, while minimizing oxidation due 
to trapped air. 

Containment Study 

The experimental set up utilized in this study 
is illustrated in Figure 3. 

A tube furnace constructed by winding Kanthal 
resistance wire around a high grade alumina tube and packed 
with refractory firebrick insulation is used. A high 
density mullite tube is mounted within the alumina tube 
to provide a contamination-free environment (with minimum 
outgassing of the heater tube) and to provide more uniform 
temperatures in the effective reaction zone. A gas inlet 
port is cemented at one end of the mullite tube. At the 
other end a nearly gas tight closure has a sliding face 
shield to facilitate insertion and removal of the specimen 
holder and temperature measurement. Tantalum baffles are 
inserted into the mullite tube between the specimen 
holder and the ends to minimize heat loss to the ends 
and to provide a uniform temperature zone over the length 
of the specimen holder. The construction of the specimen 
holder is illustrated in Figure 4. The entire holder is 
constructed of 304 SS and is free of warpage upon heating and 
cooling. The alloy /base material samples are placed in 
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individual receptables in the holder and are removed by 
means of a tight fitting SS304 strip bent to fit the 
receptacle. Each sample is further isolated in its 
receptacle by means of a tight fitting SS304 sheet which 
slides over the holder, thus preventing excessive 
volatilization and loss of elements and contamination of 
samples. 

The reaction temperature is maintained 100-200°C 
above the melting point of the alloy, to allow enhanced 
reaction/dif fusion/penetriition rates and ensure complete 
melting of the alloys. Volatilization of the elements from 
the alloy is much lower than it is from the elements 
themselves and does not require using the minimum temperatures 
required for melting. 

The specimen holder with the samples isolated by 
the covering strip is introduced into the furnace which is 
flushed with argon and maintained at temperature. 

The furnace is closed off, and an argon flow sufficient 
to maintain a pure atmosphere is provided. The durations 
of the exposures used are 24, 48, 100, 300 hours, done in 
individual batches. Following each heat the specimens 
are cooled to room temperature in the argon atmosphere by 
switching off the furnace heat. The samples are then 
removed carefully and identified prior to analysis. 


Analytic Procedures 

The samples obtained in the containment study 
(alloy element or alloys on base material) are visually 
inspected for extent of wetting and reaction. 

In those instances that the alloy wets the 
base or forms a sessile drop on cooling, the effect 
is recorded using an optical metal lograph. 

The alloy is then ground off the base material 
to within 0.5mm of the surface, and the base is sectioned 
by means of a diamond impregnated cutoff wheel. 

Different grades of abrasive wheel are used for the SS304, 
316 and Sic samples. The cutoff position is selected 
to obtain a representative section. 

The sample is then mounted in epoxy resin (cured 
at room temperature) . Alumina fillers are used to improve 
the hardness and abrasion resistance of the mold to more 
closely match the sample, thus providing better edge 
retention and minimizing relief between sample and mount. 

The samples are then ground on Sic papers down to 
600 mesh and subsequently polished on diamond impregnated 
cloths (low nap for better retention of inclusions and 
precipitates) for SiC and SS304 and 316. Alumina polishing 
abrasives are used also in the case of the stainless steels 
The final grade of abrasive used is down to 0.25pm. 



The sample surface is thus polished flat and free 
of surface relief, well suited for SEM observation and 
EDAX analysis. 

Problems associated with smearing the softer 
alloying elements are thus minimized and practically 
eliminated by grinding off excess alloy, sectioning 
through the unattacked rear face of the sample, grinding 
and polishing while keeping the specimen reaction edge 
suitably oriented, and cleaning efficiently by means 
of swabbing with alcohol, using ultrasonic baths between 
preparation stages and at the end of preparation. 

Etching is only carried out on as-received samples 
of the base materials and on base materials subjected 
to the same thermal cycles under argon as the alloy/base 
material samples, to remove artifacts due to thermal 
cycling alone, such as alloying element segregation. 

SEM micrographs are used to study the growth of 
surface reaction layers and penetration-infiltration of the 
base material, while EDAX-electron microprobe analysis, 
performed on beam-size spots and in-line scans at successive 
depths into the base material from the reaction layer , 
are used to study elemental diffusion and distribution. 
Suitable corrections (ZAF, magnification, specimen tilt 
with respect to electron beam and detector window) are 
applied and spectral artifacts eliminated prior to analysis 
of the results obtained. 


This system of analysis is seen to provide an 
effective, reliable means of studying the containment 
problem. 


Preliminary qualitative results are included in 
Figures 5a-f. 


SUMMARY 


The alloy systems used are listed in Table 1 along 
with related physical and thermal properties. 

The base materials studied include SS304, SS316L 
and SiC ("Carbofrax" from Carborundum Co.). 

- Results obtained to date indicate considerable 
reactivity of the stainless steels towards several of the 
heat storage alloys, while the SiC remains inert or forms 
a very slowly growing reaction layer with the molten 
alloys. The surface oxide on SiC has been found to be 
principally responsible for the poor wetting characteristics 
exhibited by the alloys. However, since several examples 
are cited in the literature wherein the liquid metal 
finally diffuses in through a resistant surface layer to 
attack the bulk or undermines the surface layer by initial 
infiltration at a crack, additional studies are being 
undertaken wherein the SiC^ layer on the SiC samples is 
removed to the extent possible by means of treatments with HF. 


Figure 5. Reaction of stainless steels and silicon carbide 
with molten heat storage alloys. Elemental 
diffusion curves are based on the longest 
exposure indicated. Concentrations are relative. 
The apparent effective thickness of the reaction 
layer and its dependence on time is shown at the 
top of each curve. 
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Additional areas of research worthy of study in 
finding a solution to the present problem include: 

a) Special surface or heat treatment of the 
stainless steels. 

b) Studies on SiC surface layer structures and 
the dependence of liquid metal interaction on this aspect. 

c) The development of SiC coatings and their use 
with stainless steels or other suitable structural material 
to overcome certain commercial drawbacks in the shaping 

of complex SiC shapes. 

d) Containment studies carried out over long 
periods and under pressures such as that to be found at 
the bottom of a tall liquid alloy bath. 

e) Effect on mechanical properties of the base 
materials due to reaction/infiltration by liquid metals. 
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